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Abstract 

We examine the distance of the two galactic microquasars GRO J 1655-40 and A0620-00 which are potentially the 
two closest black holes to the Sun. We aim to provide a picture as wide and complete as possible of the problem of 
measuring the distance of microquasars in our Galaxy. The purpose of this work is to fairly and critically review in 
great detail every distance method used for these two microquasars in order to show that the distances of probably 
all microquasars in our galaxy are much more uncertain than currently admitted. Moreover, we show that many 
confirmations of quantitative results are often entangled and rely on very uncertain measurements. We also present a 
new determination of the maximum distance of GRO J 1655-40 using red clump giant stars, and show that it confirms 
our earlier result of a distance less than 2 kpc instead of 3.2 kpc. Since it then becomes more likely that GRO J1655^40 
could originate from the stellar cluster NGC 6242, located at 1 .0 kpc, we review the distance estimations of A0620- 
00, which is so far the closest black hole with an average distance of about 1 .0 kpc. We show that the distance methods 
used for A0620-00 are also problematic. Finally, we present a new analysis of spectroscopic and astrometric archival 
data on this microquasar, and apply the maximum-distance method of Foellmi et al. (2006). It appears that A0620-00 
could indeed be even closer to the Sun than currently estimated, and consequently would be the closest known black 
hole to the Sun. 



1. Introduction 

Microquasars are stellar binaries in which one of the 
component is a black hole. The stellar companion is fill- 
ing its Roche lobe, ejecting matter through the first La- 
grangian point, which then organizes itself around the 
compact object as an accretion disk. The temperature of 
the disk increases strongly towards the center and large 
amounts of X-rays with powerful, persistent and colli- 
mated jets are produced. In some cases, the relativistic 
jets appear to be superluminal. Microquasars are galac- 
tic laboratories of high-energy phenomena, and they 
must be seen as part of a large paradigm where AGNs, 
microquasars and possibly gamma-ray bursts (GRBs) 
share similar physics ( jMirabel 20041. 

There are various physical reasons why the distance 
is an important parameter in our understanding of mi- 
croquasars. Firstly, the superluminal velocity effect of 



the jets (which is a relativistic time delay effect pro- 
ducing an apparent velocity on the sky larger than the 
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speed of light; see e.g. Rees 1966) depends on the dis- 
tance. Although it is a pure geometrical effect, it has 
a minimum threshold of f3 = v/c > 2~^ 2 (with c the 
speed of light) below which the jets are not superlumi- 
nal, whatever their true spatial speed v. When one tries 
to model the jets, and more importantly, the jet's launch- 
ing mechanism, it is important to know the true velocity 
(and hence power) of the jets. Moreover, this effect is 
not very common, and it will become more important to 
correctly identify superluminal objects when the statis- 
tics will grow with future X-ray missions. 

Secondly, the possible misalignment of the jets with 
respect to the accretion disk axis. As a matter of fact, 
GRO J 1655^-0 is often quoted as the typical micro- 
quasar where the jets and the disk are misaligned (e.g. 

2002| l. This is directly related to the fact 



Maccarone 



that the determination of the jet projection angle is 
linked to that of the distance, as explained below. The 
possible misalignment of the jet is an important clue of 
the formation and evolution of the system (see e.g. |Mar-| 
|tin etal.||2008] l. 

Thirdly, the black holes' origin, and their orbit in 
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our Galaxy. This is particularly true for GRO J 1655- 
40, where a change in the distance by a factor of three 
makes a large difference to the orbit of the object in the 



Galaxy. See for instance Mirabel et al. ( 2002 1 who have 
calculated the orbit for GRO J 1655-40 with D = 0.9 
and 3.2 kpc. This is important for our understanding of 
the origin of black holes in our Galaxy, and their distri- 
bution throughout the disk. 

Finally, the true luminosity of the stellar companion 
can also be an important point in our understanding of 
the microquasars as dynamical objects, and in particu- 
lar how the companion is affected by the filling of its 
Roche lobe and is certainly (partially) irradiated by X- 
rays from the accretion disk (e.g. Dub us etaL] [T999). 

The distance of microquasars is a simple yet central 
parameter in our description of these objects, condition- 
ing a large part of our understanding both of the physics 
and the astronomical views of black hole stellar sys- 
tems. 

1.1. The origin of the disagreement on the distance of 
GRO J1655-40 



The origin of this work is to be found in the publi- 
cation by Mirabel et al. ( 2002} who note that the dis- 
tance of GRO J1655-40 can be actually radically dif- 
ferent from the accepted distance of 3.2 kpc determined 
by |Hjellming and Rupen| ( |1995| l. This later value has 
since then been apparently confirmed by many other 
studies. In fact, GRO J 1655^-0 is a runaway black hole 
and Mirabel et al. published its proper motion, obtained 
with the Hubble Space Telescope. The opposite direc- 
tion of the proper motion points almost perfectly (see 
their Fig. 1) to the center of a cluster (NGC 6242) lo- 
cated at 1 .0 kpc from the Sun ( Kharchenko et al. 2005 ). 



It is then tempting to think that the system received a 
kick velocity at the moment of the primary's supernova 
explosion, and then moved away from its original clus- 
ter. This is the starting question: is the distance of 
GRO J1655^0 3.2 or 1.0 kpc? Among other things, 
the jets are superluminal with D — 3.2 kpc, but not at 
1.0 kpc. 



Foellmi et al. (2006) published a new method pro- 



viding only an upper limit to the distance and giving 
D < 1.7 kpc when applied to GRO J 1655^-0, thus chal- 
lenging the accepted distance of 3.2 kpc. This method 
is also partially problematic, and it will be discussed 
critically below. With a maximum distance of 1.7 kpc, 
GRO J1655^-0 becomes a likely candidate of being the 
closest (stellar) black hole to the Sun. The position is 
currently held by A0620-00, which has an average mea- 
sured distance of about 1.0 kpc. We will show however 



that the distance of A0620-00 is also very uncertain, re- 
vealing different problems in each distance method. 

In summary, none of the published distances of 
GRO J1655^10 and A0620-00 are so far reliable. This 
paper emphasizes the difficulties encountered when de- 
termining the distance of microquasars, and shows that 
the distances of probably many, if not all, microquasars 
in our Galaxy are much more uncertain than currently 
admitted. In this paper we aim to provide a fair and crit- 
ical review on the distance methods used on these two 
galactic microquasars. Other microquasars will be the 
subject of a subsequent work. In addition to |Foellmi 



et al. (2006), partial material (incomplete and par- 



tially incorrect on the radio-jet distance of GRO J1655- 
40) has already been published in [Foellmi| ( |2006[ ) and 
Foellmi { 2001) . This paper addresses many more im- 



portant details, and in particular the issues of the dynam- 
ical studies used to determine the absolute magnitude 
of the secondary star of GRO J 1655^-0. Moreover, we 
also provide a completely new estimate of the distance 
of GRO J 1655-40. Finally, we apply the method of 



Foellmi et aL] ( |2006) to A0620-00, along with archival 
astrometric data. We conclude that GRO J 1655^-0 is 
certainly closer than the current accepted distance, but 
that A0620-00 might be even closer still. 

1.2. Distance methods: the basics 

The distance of a stellar object is often measured by 
comparing its absolute and apparent magnitudes. There 
are other methods, such as astrometry and parallax, and 
jet speed measurements (see below) for the special case 
of microquasars. But most often, the core method is 
simply this one. It reads: 



- M = 5 log(D) - 5 



(1) 



where m aae is the true apparent magnitude of the object, 
affected only by the geometrical distance separating the 
object and the observer. 

There are (at least) two major issues with this method: 
one is general, and one is specific to stellar binaries of 
short period. The first observational difficulty is that 
m tme is usually different from the observed apparent 
magnitude m ^ s since the light is going through some 
absorbing patchy interstellar medium, which makes the 
star dimmer and redder. We have: m true = ;« b s - A, 
where A is the absorption, in the given passband, ex- 
pressed in magnitudes. Therefore: 



m ohs - M - A = 5 log(D) - 5 



(2) 



The determination of A is critical. It is often not mea- 
sured directly, but rather is the so-called color excess, 



or reddening: E(B - V), which is the relative amount 
of additional red color between the B and V bands due 
to the absorption of bluer wavelengths by the gas and 
the dust. For a star, the color excess can be obtained by 
comparing the observed and the intrinsic B and V color 
indices (the latter being inferred from the spectral type): 



E(B - V) = (B - V) - (B - V) 
The absorption in the V band follows: 
A v = R x E(B - V) 



(3) 



(4) 



where R is the total-to-selective absorption. For other 
colors, this relation needs corrections (see e.g. |Fitz 
patrick) |1999| l. The practice shows that what is often 
measured directly is E(B - V), R being often approx- 
imate^ by a value between 3.0 and 3.7 with a more- 
or-less canonical value of 3.1. In fact, R is a function 
of the interstellar reddening curve and the color of the 
stars, because the wide passbands of the photometric B 
and V filters makes the effective wavelengths of the fil- 



ters shift with different stellar intensity gradients. Olson 



( 1975 i gives an approximate relation for R: 



R = 3.25 + 0.25 (B - V)o + 0.05 E(B - V) (5) 

where (B-V)o is the unreddened color index of the star. 
This relation is valid for normal stars with (B-V)o < 0.4 
and E(B - V) < 1 .5 to within an error of 0.05 in R dOT 



son 



T975J. See for instance Crawfor d and Mandwe 



wala (19761 for a comparison for various photometric 
systems, and McCall (2004) for an updated discussion 
on this topic. 

The other issue with the method consisting of com- 
paring the magnitudes, specific to the stellar binaries 
with short period, is that the stellar companion (here 
the secondary star) is certainly not spherical anymore, 
since it completely fills its Roche lobe to feed the accre- 
tion disk. Moreover, its surface temperature might also 
not be homogeneous because of irradiation (e.g. Dubus 
et al. 1999| O'Brien et al. 2002), making its average 
observed temperature a function of the orbital phase. 
Hence uncertainties arise when one tries to estimate the 
true value of the absolute magnitude, M, in equation [2] 
since it is often calibrated with "normal" (spherical) and 
isolated stars. 

The absolute magnitude is normally obtained through 
the determination of the spectral type of the star, which 
gives the temperature r e ff, and the modelling of either 
the radial-velocity curve or the multi-color lightcurves 



'The difference between R =3.0 and 3.7 implies a relative distance 
uncertainty of 1.38. 



of the binary system gives the size of its orbit. From the 
latter, one can compute the effective Roche-lobe radius 
(see for instance Eggleton[ |1983[ ), which is then identi- 
fied to the radius of the star. Assuming a uniform tem- 
perature distribution across the star and that it is roughly 
spherical, the luminosity (and thus the absolute magni- 
tude) can be estimated with: 

L = 4na-R 2 T 4 (6) 

where <j is the Stefan-Boltzmann constant. 



2. The published distance of GRO J1655-40 

GRO J1655-40 (a.k.a. Nova Sco 94) has been dis- 
covered as a Soft X-ray Transient (SXT) on July 27, 
1994 with BATSE on board the Compton Gamma Ray 



Observatory ( Zhang et al. 1994 1. Its jets were observed 
in the radio, giving a distance of about 3 kpc, imply- 



ing that they are superluminal: 1.5±0.4 c (Tingay et al 



1995), 1.05 c (Hjellming and Rupen 1995 1, where c 



is the speed of light. GRO J 1655^-0 was the second 
superluminal source in our Galaxy shortly after the dis- 



covery of the jets in GRS 1915+105 by Mirabel and 



|Rodriguez|(|1994[). The value of 3.2 kpc has been pub 
lished by Hjellming and Rupen (1995) based on VLA 
and VLBA radio data. 

Since then, many publications use directly this canon- 
ical distance. For instance: Brandt ~et al.| ( 1995) ; Barret 
|et al.|(fl99"6l);|Reg6s et al.| (|1998)>; |van der Hooft et al. 
p9^pmlips et al.l(|1999f;|Shahbaz et al.|(|1999|);|Kvr 



lulkers et al.| ( |2000| i; |Soria et al.| ( |2000| i; |Greene etaL 
( |27X)!) ; |Combi et al.| ( [200TT ); |Buxton and Vennes| ( |200TT ); 
|Yamaoka et al.| ( |2001[ ); |Gierlirisk i~et al.| ( |2001[ ); |Kub 



ota et al. (2001 ); O'Brien et al. (2002); Remillard et al. 



(|2002l>; |Kong et al.| l |2002); |Kobayashi et al.| ([2003]>; 
Stevens et al.| | |2003) ; |Willems et al.| fl2005]l; |Brock-' 



soppetal.|( 2006); Miller et al.|(|2lX)6) ; ICaballero-Garcia 
et al.| ( |2007| l; |Martin et al.| ( |2008| l and fChakrabarti et~ 



(2008) . More recently Shaposhnikov and Titarchuk 

(2009) use GRO J 1655^-0 and its parameters as a ref- 



erence object for determining the mass of Cyg X-l. 
But a significant number of studies seem to confirm 
the canonical distance. We challenge here not only the 
lower limit of 3.0 kpc of GRO J 1655^40 but also every 
published confirmation. 



We note however that Migliari et al. (2007) use with 



caution the upper limit of 1 .7 kpc determined in Foellmi 
[eTaLl ( |2006l l. 

2.1. The radio-jet kinematic distance of 3.2 kpc 

|Hjellming and Rupen| (fl995 ) present new radio data 
from which they infer a value for the distance of 



GRO J1655-40. The method is simple: the opposite 
motions of the receding and approaching jets and fi + 
respectively) are directly related to jet projection angle 
relative to the line of sight 9, the true jet speed B = v/c, 
and the distance D. The authors measure the intrinsic 
proper motions of each jet: 54 mas d -1 for the NE com- 
ponent, and 45 mas d _1 for the SW component, and use 
the kinematic equation described in Mirabel and Ro- 
driguez] ([1994): 



Hi 



/?sin(0) 



1 + Bcos(9)D 



(7) 



where c is the speed of light. There are two equations 
and three unknowns. A constraint on the distance can 
be obtained by eliminating 9 and requiring that B < 1 : 



D I 2/i+jU- 



C \/J.+ +jU- 



< 1 



(8) 



Taking v - c, the maximum distance of GRO J 1655- 
40 inferred from the proper motion of the radio jets is 
D < 3.5 kpc. 

A constraint on the inclination angle of the jets can 
also be obtained by eliminating v/D and requesting that 
V < c. We obtain: 



! = OfzZifi) cos(0) -i <1 

C (yU_+// + ) 



(9) 



This gives: 9 < 84.8°. Rearranging equation|7]and elim- 
inating B, we can write, using the value of 9: 



1 



2p + p.- 
v/D + p_ 



= sin0 < 0.996 



(10) 



that is v/D > 49.3 mas d~'. All these results are given 
in |Hjellming and Rupen| ( |1995[ ). 

At this point the authors mention that: "For a dis- 
tance of 3.2 kpc, this corresponds to V > 0.91c, implying 
84.3° < 9 < 84.8°." Why 3.2 kpc? The only previous 
mention of the distance at the beginning of the article is 
only stating that the source lies "at a distance of about 3 
kpc", and for which three references were given (all dis- 
cussed below ): [ McKay and Kesteven ( 1 994) > ; Harmon 
et al.| ( |1995| >; |Tingay et al.| ( |1995) . In fact, the value 
D = 3.2 kpc is constrained to lie between the value 
of 3.0 kpc considered as a lower limit and determined 
elsewhere, and the actual result on the upper limit of 3.5 
kpc. However, choosing 3.2 kpc fixes 9 that is now con- 
strained to a very narrow range (0.5°), and the system is 
now defined. 



2.2. The confirmation of 3.2 kpc with possible wiggles 
inferred from unaligned images 
Furthermore, this value of 3.2 kpc is being confirmed 
by a model of the possible jet precession (adding two 
parameters to the kinematic model: the precession pe- 
riod and axis inclination). It is however barely relevant 
concerning the distance, since, as mentioned earlier, the 
systems characteristics are now fixed. This additional 
modelling gives 9 = 85° based on the fact that the jet 
precesses, or more precisely: "[...] the jets 'wiggle' 
slightly about the best-fit position angle". The authors 
performed a detailed modelling, similar to what has 
been done by one of the authors on the galactic source 



SS433 which shows clearly jet precession (Hjellming 
land Johnston| 1 1 98 1 a|b[ [198^ . 

The model of the wiggles assume, of course, that they 
are true, and are due to the kinematics only. In Fig. 3. 



of Hjellming and Rupen (T995) we see the variations 
of position of roughly 4 different ejectas (following the 
main solid lines only). It is obvious that few points per- 
fectly follow the best-fit constant proper motion. But 
it is not obvious at all that they form a regular periodic 
pattern. Succeeding at modelling these wiggles does not 
mean we can interpret them as the signature of jet pre- 
cession taking place in GRO J 1655^-0, since the low 
number of points makes the fit poorly constrained. At 
that point, it should be clear enough that these supposed 
wiggles do not help to secure the distance. But there are 
additional problems attached to them. 

In order to interpret wiggles as precession, one must 
ensure that the global motion follow lines of constant 
proper motion. But the 22 epochs of VIA observations 
used for by this model did not resolve the source at a 
level of 100 mas, but only as a multi-core object elon- 
gating with time. The reason why a constant proper mo- 
tion is a reasonable hypothesis is because it is consistent 
with what is seen in the VLBA observations. But the 
authors emphasize the lack of a very-long-baseline in- 
terferometry calibrator, which implies that these VLBA 
data are self-calibrated, "eliminating all absolute posi- 
tional information, and leaving the alignment of the dif- 
ferent images a free parameter." Consequently, the fact 
that the brightest point in each image is the stationary 
center of ejection is an hypothesis. 

Moreover, as mentioned in the paper, "these [VLA] 
and other [unspecified] data are consistent with constant 
intrinsic proper motion of 54 mas d _1 " and later "the un- 
derlying proper motions appear constant". This value is 
in agreement with the result of equation 10 The hy- 



pothesis of constant proper motion is strengthened by 
the fact that the daily Southern Hemisphere VLBI Ex- 



periment (SHEVE) array observations of Tingay et al. 



(1995) are consistent with the major structures of the 
VLBA images of Hjellming & Rupen. As a matter of 
fact, the authors note that the proper motion inferred 
from SHEVE data of 65 ±5 mas d -1 actually agrees 
with the 62 mas d -1 motion of the outer edge of the 
early NE ejecta. All these measurements indeed appear 
to roughly agree, but the uncertainties are likely large 
enough to encompass the wiggles. In summary, not only 
are the mere existence of these wiggles doubtful, but the 
simple ability to extract meaningful and additional con- 
straints from it is also questionable. 

We conclude that the value of 3.2 kpc has been cho- 
sen between a firm upper limit of 3.5 kpc and the exter- 
nal indication that its lower limit must be "about 3 kpc", 
which consequently fix the value of 6. A questionable 
model is used to strengthen its value, and, as an obvious 
consequence, confirm the distance value. Literally, later 
in the text, it is said that the "[...] kinematic model for 
GRO J 1655-40 gives a distance of 3.2+0.2 kpc". 

2.3. Where does 3 kpc come from? 

As noted above, three references are given for a first 
estimation of the distance: Harm orTet al.| ( |1995| l, |Tin-| 
|gay et al.| ( fl995| l and |McKay and Kesteven] ( |1994| i. Un- 
fortunately, the first reference is literally citing the two 
others for the distance value, and must therefore be dis- 



carded. The paper by McKay and Kesteven ( 1994! is 



actually an IAU Circular which simply states that "HI 
observations of GRO J 1655-40 made with the AT Com- 
pact Array show solid absorption in the velocity range 
+ 10 to -30 km s _1 , with a further isolated weak feature 
at —50 km s . The balance of probabilities is that the 
distance is around 3.5 kpc, unless the —50 km s fea- 
ture is due to an atypical cloud." Although not being a 
robust measurement, the result obtained in this Circular 
needs to be verified. 

Identifying the origin of the negative-velocity fea- 
tures in the absorption spectrum is the crucial point, 
since it provides an estimate of the lower limit of the dis- 
tance, if correctly interpreted (i.e. if correctly identified 
and attributed to a component whose velocity can be es- 
timated). To interpret the radio spectrum, one must con- 
sider a background source of unknown distance emit- 
ting continuum radiation (here GRO J 1655^-0) that is 
being intercepted by foreground HI clouds. When one 
looks inside the solar galactic orbit, the line of sight 
goes through the multiple galactic spiral arms. Assum- 
ing all the clouds are moving with the mean galactic ro- 
tation scheme, the more distant from the Sun the cloud 
is, the more negative is its velocity, in the Local Stan- 
dard of Rest, down to the tangential point where the dis- 
tance/velocity relation flips back. 



Tingay et al. ( 1995 1 presented new radio VLBI and 
ATCA data of GRO J1655^10. It is the only true work 
studying the lower distance limit of GRO J 1655^-0. 
Firstly, we note that a rather large distance is expected 
by the authors in order to agree with: "a significant red- 
dening due to absorption' 
p994) i. |della Valle| ( [1994) is a IAU Circular stating 
nothing more than: "[...] The [optical] spectrum ex- 
hibits prominent, broad Balmer lines [...] superimposed 
on a relatively red continuum. [...]" However, the spec- 
trum has been taken during a flaring state. The spectral 
type of the secondary, and its possible veiling by the ac- 
cretion disk, were, at that time, unknown. Leaving this 
aside, let us concentrate on the radio spectrum. 

The HI spectrum of Tingay and collaborators, ob- 
tained with ATCA (see their Fig. 2), shows a multi- 
component profile, with strong features at — h5 km s _1 , 
between -10 and -20 km s _1 , and isolated weak fea- 
tures at -30 and -50 km s~' (-18% and -2% of the 
normalized continuum flux respectively). It is said that 
the latter feature is confirmed with multiple observa- 
tions but no references are given. 

According to Tingay et al. ( 1995| l, the feature at 
-50 km s _I would imply a lower limit of ~4.2 kpc if 
it was participating to the mean galactic rotatiorj^] This 
feature is actually discarded by Tingay and collabora- 
tors because it cannot be ruled out that such feature is 
driven by an expanding shell surrounding the Scorpius 
OBI association located at 1.9 kpc from the Sun (see 
Fig.|TJ. This is however the feature used by McKay and 



Kesteven ( 1994 1 to derive an approximate distance of 



3.5 kpc. In other words, McKay and Kesteven (1994) 
derive D ~ 3.5 kpc thanks to the feature at -50 km s" 1 
that is discarded by Tingay et al. (1995]) because it 
would imply D > 4.2 kpc. 

According to Tingay et al. (1995]), the feature at 
-30 km s -1 implies a lower limit of 3.0 kpc if it is as- 
sociated with the mean galactic rotation. To strengthen 
their conclusion, they compare their spectrum with that 
of a nearby region GRS 345.4+1.4 (a.k.a. CTB 35 A) 
studied by Caswell et al. (1975]) and located 2.4 kpc 
away. In fact, the lower limit of the distance of 
GRO J 1655^1-0 is built on this hypothesis: the feature 
at -30 km s _1 in the radio spectrum is moving with the 
mean galactic rotation scheme, which is confirmed by 
the comparison with GRS 345.4+1.4. 

But if there is such uncertainty on the interpreta- 
tion of the feature at -50 km s _1 , why can't the 



2 We assumed that the authors follow the same galactic model of 
|Caswell et al.|ft975) which is not explicitly stated but whose results 
are explicitly used. 
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Figure 1: Map of the region around GRO J 1655—40, where the re- 
gion GRS345+1.4 studied by Caswell et al. ( 19 75) is also shown as 
a big open circle. The 24 stars studied by Crawford et al. 1 19891 are 
indicated (stars belonging to Sco OBI in blue circles, field stars in 
black squares with names). It shows that GRO J1655^10 is located 
at the core of the region studied by Crawford and collaborators. Note 
the presence of the star HD 152269 which lies close to the center of 
Sco OB 1 but is a foreground star. The green triangle to the left shows 
the position of the comparison F6IV star HD 156098 used by Foellmi 
|et al.| j2006 1 in their analysis. 



feature at -30 km s" 1 not also be associated with 
Sco OBI? Or could the latter be indeed associated with 
Sco OB 1 while the most negative one is not? The spec- 
trum of GRS 345.4+1.4 looks indeed similar to that of 
GRO J 1655^-0, except that it has no absorption feature 
with velocities more negative than -24 km s _1 . But as 
shown on the map in Fig. [T] GRS345.4+1.4 is closer 
to Sco OBI than to GRO J 1655^10. Why does the 
spectrum of GRS 345. 4+ 1.4 shows nothing more neg- 
ative than -24 km s _1 although it is angularly closer 
but still behind the association compared to the position 
of GRO J 1655^0? 

( fT975l > confirm 



We note moreover that Caswell et al. 



the results of Radhakrishnan et al. ( 1972) that the HI ab- 



sorption spectrum of the region GRS 345.4+1.4 is prob- 
lematic since either the foreground absorbing cloud at 
V = -24 km s~'or the HII region behind it (which is re- 
sponsible for the continuum emiss ion) ha s a "pe culiar" 
motion of ~10 



km 



Caswell et al. (1975) deter- 



mine a (near kinematic) distance of 2.4 kpc by assum- 
ing a mean velocity of the multi-profile HI absorption of 
~ -20 km s _1 . An uncertainty of 10 km s^'translates to 
an uncertai nty of about 1 .5 kp c in distance, if on e uses 
the scale in |Caswell et al. (19750 Furthermore , | Shaver) 



|et al.| ( |1982| l have shown that there are a number of HI 
clouds in this region that have a peculiar velocity which 
cannot be accounted for by assuming that such cloud is 
moving with the mean galactic rotation, since the de- 
rived distance from the HI absorption appears larger 
than that derived from optical observations of HII re- 
gions along the same line of sight. We are therefore 
entitled to conclude that no meaningful comparison can 
be made between two absorption radio spectra in this 
region of the sky and the interpretation of the negative 
velocity features in the absorption radio spectrum is un- 
certain. 

We conclude that the lower limit of 3.0 kpc on the 
distance of GRO J 1655-40 has not been established. 

2.4. Optical spectroscopy of stars in the direction of 
GRO J1655-40. 



|Crawford et al.| ([1989) have studied the interstellar 
sodium and calcium absorption lines towards the Scor- 
pius OB 1 association with high-resolution optical spec- 
troscopy. The stars observed in this study are also 
shown in Fig. [T] They observed that all the spectra 
of Sco OBI members show structured features in the 
sodium doublet lines with velocities spanning a range 
as wide as 40 to 60 km s , to the contrary of all other 
field stars, including the foreground star HD 152269 lo- 
cated right at the center of Sco OBI on the sky. The 
absorption features at the bluer wavelengths inside the 
absorption lines are interpreted as truly blueshifted be- 
cause the spectra are calibrated with atmospheric water 
lines directly on the spectra (i.e. the zero-point in veloc- 
ity is known). 



Crawford et al. ( 1989 1 made the global follow- 



ing conclusions, despite variations between individual 
sources. Firstly, absorption lines comprised in the range 
-20 < v'heiio ^ km s must arise from material be- 
tween Sco OB 1 and the observer, since it cannot partic- 
ipate to the mean rotation curve. Second, the absence of 
blueshifted absorption features with Vh e i < -20 km s _1 
in the spectrum of HD 152269 implies that the absorb- 
ing material responsible for these blue features is com- 
prised between the star (720 pc away), and the clus- 
ter. Given an average velocity of Sco OBI of about 



3 We note that Caswell and coworkers use actually the galacticj 



model of Schmidt 1 1965 1 with a galactic center distance of 10 kpc, 
while the modern value is about 8.0 kpc (Groenewegen et al.||2008} . 
The extraction of the velocity-distance relationship directly on the fig- 
ure of Caswell et al's paper reveals that the printed scale does not 
correspond exactly to the results written in the text itself (a velocity of 
-24 km s~'correponds to 2.8 kpc instead of 2.4 kpc; possibly causing 
the confusion). However, the use of a smaller galactic center distance 
roughly compensate this error on the velocity-distance relationship. 
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Figure 2: The sodium doublet observed in the spectrum of 
GRO J1655^K) as function of the heliocentric velocity. Rest wave- 
length of the sodium were used: /U5889.951, 5895.924. The spectra 
were taken from Foellmi et al. ( 2006}, and were combined together 
without preliminary radial-velocity shift. The zero-point of the ve- 
locity scale in the figure is simply that of the rest wavelength of the 
sodium lines, to which the calibration of the spectrum agree within 
5 km s (see text). Four gaussians were fitted to each independent 
line. The top panels show the spectra (points), with the individual 
gaussians in color, and the total gaussian in gray line. The bottom 
panels shows the difference between the spectrum and the fit. 



-25 km s" 1 , the authors conclude that the absorptions 
seen at the bluest wavelengths of the sodium doublets 
are caused by expanding material associated with the 
cluster. 

In order to compare with GRO J 1655^-0, we have 



taken the original spectra published in |Foellmi et al. 



( 2006 ), and averaged them with no radial- velocity shifts 
prior to the combination, in order to build a mean spec- 
trum where the (static) interstellar lines are well aver- 
aged. The spectrum is shown in Fig. [2] Our UVES spec- 
tra are well aligned in a relative manner, since all sharp 
reddest wings of the Na lines appeared well aligned 
to each other. Moreover, we have identified about 10 
(residual) sky lines inside the averaged spectrum (also 
identified in Foellmi et al. 2006). When compared to 



the skylines atlas of |Osterbrock et aL (1996i, the ve- 
locity zero-point our the spectrum appears to be correct 
within 5 km s _1 . We have fitted the interstellar lines 
with four gaussians. We note the differences, in veloc- 
ity, between the bluest and the reddest gaussians are 5 1 
and 52 km s respectively for the two lines. 

We can see in Fig. [2] that there is no absorption be- 
tween and ~ -30 km s _1 . Along with the conclu- 
sions of Crawf ord et al.| ( fl989} , it would mean that there 
is a very small amount of absorbing material between 
GRO J 1655^-0 and the observer. This conclusion is 



hardly reconcilable with a large distance (and afortiori 
larger than that of Sco OBI). Moreover, the sodium 
absorption is saturated between -20 and -60 km s _1 
and a weaker blue feature exist at ~ -85 km s . It 
seems more likely that these absorptions arise from an 
expanding material related to GRO J 1655^-0 itself. As 



matter of fact, Combi et al. (20071 have investigated 
the region of GRO J1655^W with radio and infrared 
data tracing the HI, CO, and the gas morphology of 
the region. They have found that there is evidences 
of an HI hole of l°.5xl°.5 in diameter at a distance of 
D — 1.2 + 0.4 kpc, and compressed CO material accu- 
mulated in part of the shell border, as well as infra-red 
emission with characteristics of shocked-heated dust. 
This study strengthen the suggestion that GRO J 1655- 
40 could be located around 1.0 kpc, and possibly origi- 
nated from NGC 6242. 

3. The difficulties in estimating the extinction to- 
wards GRO J1655^10 

After the value of 3.2 kpc had been published, a 
significant number of studies confirmed or strength- 
ened this result, by using various distance methods. 
We discuss in the following sections the various prob- 
lems encountered when estimating the extinction to- 
wards GRO J 1655-40. 

3.1. The validity range of the sodium equivalent width- 
color excess relationship 

A relationship exists between the equivalent width 
of the sodium doublet or the calcium lines in an op- 
tical spectrum, and its color excess. This relation- 
ship has been used, for instance, by Bianchini et al. 



( 1997} for GRO J 1655-40, although using spectroscopy 
of low resolution preventing them to see the saturated 



and multi-profile nature of the lines (see Foellmi et al 



|2006| for the details, and Fig.|2]above). They use the re- 
lationships given by Herbig ( 1975 ) for the c alcium line 
at 6613A (see Fig. 4 of Herbig 1975) and |della Valle 
and Duerbeck ( 1993} (wh o actually use the photoelec 



trie photometry of Cohen 1975 1 for the sodium doublet 
at 5980A: 



E(B -V)~ 0.61 x £W NaI _ D - 0.08 (11) 
Bianchini et al. ( 1997} obtained a color excess rang- 



ing from 0.97 to 1.30 mag, from which they adopt a 
mean value of 1 . 1 3, and state that it is in agreement with 
the value found by other studies on this target (for e.g. 
|Bailyn et al.| [1995a but see below). But such a range 
on the color excess translates, using a standard value of 



R — 3.1, to an uncertainty of one magnitude on the ab- 
sorption Ay and therefore a factor of 1.6 in the distance. 

Moreover, |Munari and Zwitter| (jT997j have shown 
that the unambiguous range between equivalent widths 
of Nal lines and E(B - V) is < E(B - V) < 0.4, i.e. 
much lower than the value measured by Bianchini and 
coworkers. If one used the (scattered) relation of |della| 
Valle and Duerbeck ( 1993| l, the limiting range of Mu- 
nari & Zwitter translates into equivalent widths between 
0.13 and 0.79 A i.e. with an upper limit much lower 
than the value of 2.26 A from Bianchini et al. Moreover, 
assuming a single gaussian profile, it implie^ja FWHM 
of the lines between 0.12 and 0.75 A, or a correspond- 
ing resolving power of about 50 000 for the narrowest 
Na lines, and 8000 for the broadest. It means that Bian- 
chini and coworkers simply do not have the resolution 
necessary to derive an accurate value of the color excess 
to confirm the value of 3.2 kpc. 

3.2. Six extragalactic supernovae 

A distance "of ~3.0 kpc" was proposed on the ba- 



sis of optical data by Bailyn et al. ( 1995a I. The authors 



have also measured the equivalent widths of Nal-D lines 
in their spectrum (EW = 4. 5 A) that ha s a resolution of 



10A (which is even worse than that of Bianchini et al. 



1997 and moreover blended with Hel emission). They 
finally compute a color excess E(B - V) = 1.15, us- 
ing the relation between the equivalent width and the 
color excess given by Barbo rTet al.| ( [T9 90). The latter 
paper studies the type-la supernova SN1989B in NGC 



3627. Barbon et al. ( 1990 1 determined an empirical and 



roughly linear relation between the equivalent width of 
the Nal-D lines and the color excess E(B - V) with the 
spectra of (only) six extragalactic supernovae. Their re- 
lation reads (see the end of their section §3.2): 



E(B -V)~ 0.25 x £W NaI _ D 



(12) 



defined for E(B - V) between 0.1 and 1.0, after remov- 
ing the galactic contribution. Leaving aside the intrinsic 
difficulty of determining the galactic contribution, the 
source of the data used to determine this relationship is 
not given in Barbon et al. ( 1990[ ), and we can question 



the reliability of a relation calibrated with 6 points only. 



Moreover, using the relation given by della Valle and 



Duerbeck ( 1993 i, equation 11 would give E(B - V) — 
2.66, instead of 1.15. It implies an absorption value of 
A^ pt ~ 8.24 instead of 3.56, and therefore a factor of 8.7 
in the relative value of the distance (equation El. Even 



ignoring the line saturation in the case of GRO J 1655- 
40, the two relationships are still hardly reconcilable. 

Interestingly, Bail yn et al.| ( | 1995a] > claimed that their 
result "is consistent with the EW of other interstellar 
lines in the optical domain". Finally, they use the clas- 



sical relationships of Allen ( 1973| discussed below) and 



Herbig (1975 who actually did not study the sodium 



doublet), to conclude that the distance of the source is 
compatible with D ~ 3 kpc, "in agreement with the ra- 



dio observations" of Tingay et al. ( 1995 1. 



3.3. Where is the HST/STIS spectrum? 

|Orosz and Bairyn| (fl997) presented an extended spec- 
troscopic and photometric dataset. In particular, they 
have found very clear ellipsoidal variations in their 
BVRI lightcurves, obtained in February and March 
1996, when the system was not yet completely in qui- 
escence. Although they mention various consistency 
checks throughout the paper, they do not measure the 



distance, but rather rely on that of Hjellming and Ru 



pen] ( |1995| l, said, along with E(B - V), to be "tightly 



constrained". As for the color excess, they assumed a 
value of E(B - V) = 1.3 ± 0.1, actually obtained by 
|Horne et al.| ([1996 ) who used high-quality UV spectra 
obtained with the Hubble Space Telescope. However, 
Home et al.| ([1996} is an IAU circular where the spec- 



trum is not visible, and where it is simply stated that 
"deep 220-nm absorption in the HST spectrum suggests 

E(B - V) = 1.3 mag." 

Similarly, van der Hooft et al. ( 1998 ) presented new 
VRi photometric data acquired during 28 consecutive 
nights in March 1996 with the Dutch 0.91m telescope 
in La Silla (Chile), when the source was said to be close 
to its quiescence brightness (but see below). From the 
modelling of the lightcurve they obtain a inclination an- 
gle of about 67° and a black hole mass between 6.3 
and 7.6 M Q (consistent with Orosz and Bailyn 1997| l. 
However, they have assumed again the radio distance of 



3.2 kpc by Hjellming and Rupen ( |1995| >, and a color ex- 
cess of E(B - V) = 1.3 mag, taken from |Horne et al.| 
( 1996[ l. Interestingly, they mention that the distance is 
well constrained since it has been found consistent with 



the results of McKay and Kesteven (1994), Tingay et al. 
( [T995) , |BaiTyn et al.| ( |1995a| l and |Greiner et al.| ( |1995] r 



4 FWHM Ga 



\/2 In 2 ■ o" where <x is the gaussian width. 



We have shown above how weak were these first three 
references concerning the distance of GRO J 1655— 40 
(the fourth and last reference is discussed in Sec.|4|. 

3.4. A new and lower value of the color excess of 
GRO J1655-40 
One can directly estimate the color excess of 
GRO J 1655-40. As mentioned in |Foellmi etaLlPOOS) , 



the F6IV comparison star HD 156098 has a known Hip- 
parcos distance of 50±0.2 parsec. It is close enough 
to assume that it has a negligible interstellar absorption 
(see for instance [Welsh et aL]|1990| l. Therefore we as- 
sume that its observed color index is equal to its intrinsic 
color index: (B -V) F , = (B- V) P , fi = 0.46 [using SIM- 
BAD; see also Fitzgerald ( 1970) ) who give (B - V) = 
0.46 for an F6IV]. Using the mean visual magnitude 



V = 17.12 from Orosz and Bailyn (T997JI and adopt- 
ing B ~ 18.6^from their lightcurve, the color index of 
GRO J1655-40 read s: (B - V) GRO = 1.48 mag. Since 
Foellmi et aL] ( [2006) have shown that HD 156098 repre- 
sents a fairly good twin to the spectrum of GRO J 1655- 
40^] the color excess of the microquasar simply follows: 
E(B - V) = (B - V) GRO -(B— V) Ftfi = 1-02 mag, which 
is smaller than previous values. This conclusion is also 



reached by Beer and Podsiadlowski (2002) who men 



tion that the value E(B —V) = 1 .3 is not consistent with 
a F6IV star but rather an A8 or earlier. The value of 
(B - V)f*fl is a bit outside the validity range of equa- 
tion [5] If we nevertheless use this equation, we derive 
R = 3.42 and therefore Ay = 3.49. The implications of 
these new values to the luminosity and mass ratio of the 
GRO J1655^0 are discussed in Sec.0 

4. The optical absorption derived from X-ray or op- 
tical data 

We have seen above some difficulties at determining 
the optical absorption towards a source. This optical 
absorption can in practice also be obtained from X-ray 
data. We present here the issues related to that approach 
and its application to GRO J1655^-0. In particular, we 
show that no systematic overestimation from X-ray data 
when compared to optical data can be claimed. 

4.1. Is the optical absorption of GRO J 165 5-40 from 
ROSAT flawed? 
Greiner et al. ( 1995|l presented new ROSAT X-ray 



data of GRO J 1655-40, from which they infer a dis- 
tance of 3 kpc (no uncertainties are provided). Their 
method consists of fitting the halo of the observed ra- 
dial profile of the source. This halo is produced by the 
scattering of the X-rays by the interstellar dust. To fit the 
radial profile of GRO J 1655-40 observed with ROSAT, 



There is an error in Foellmi et al. 1 12006} who use B ~ 16.65 
instead of B~ 18.6. 

6 We emphasize here that the spectra look the same, i.e. the global 
average spectroscopic parameters are comparable. We obviously do 
not mean that the stars are identical. This distinction is important in 
the next sections. 



they assume an uniform dust distribution between the 
observer and the source in their model. They obtain, 
with not many details, a value of the effective optical 
depth at 1 keV of r eff ~ 0.33. 



Furthermore, they use the results of Predehl and 
Schmitt (1995) who have studied in details X-ray ha 



los in ROSAT sources, and have shown that a good cor- 
relation exists between the simultaneous measured dust 
and hydrogen column densities: "indicating that gas and 
dust must be to a large extent cospatial". From the 
fractional halo intensity it is thus possible to derive the 



dust column density. The relation used by Predehl and 



|Schmitt| ( fl995| l reads (see the end of their Sec. 3.4): 



T el fMv = 0.056+0.01 



(13) 



where Ay is the visual absorption (expressed in magni- 
tude). 

Using this correlation, and assuming that the sight- 
line for GRO J 1655-40 has the same gas-to-dust ratio as 
the sight-lines for which the relations between Ay and 
Nh have been established, Grein eret al.| ( |199"5] l obtain 
for GRO J1655-40 an absorption of Ay = 5.6 mag (and 
a hydrogen column density of Nh = 7.0 x 10 21 cm" 2 , 
implying a color excess of E(B - V) ~1.7, assuming 
R = 3.1). This is significantly (by ~ 1.5-2 magni- 
tudes) larger than any other estimation from the studies 
in the optical (Ay ~ 3.8-4.1 mag) and our own estimate 
above. 

To finally compute the distance, Greiner et al.|(|1995 i 
use the mean extinction law given by Allen ( 1973 Sec. 
125, p. 263): Ay =1.9 mag kpc -1 , and say that it is in 
agreement with other determinations of the distance. It 
is however hard to consider the relation of Allen more 
than a basic approximation, since it is composed of two 
parts: Ay — 1.6 mag kpc -1 from "interstellar absorbing 
clouds" and Ay - 0.3 mag kpc -1 from "grains between 
the clouds". The value of 1.6 mag kpc -1 is actually ob- 
tained by simply assuming that there are 5 clouds per 
kpc in the galactic plane, and that the mean visual ab- 
sorption is about 0.3 mag per cloud (|Allen 1973 Sec. 



124, pp. 262-263). The value of 0.3 mag kpc -1 is quoted 
from Gottlieb and Upson ( 1969 ) who in fact clearly 
show, by dividing the sky into more than 200 zones, 
that the extinction in a particular direction is by far more 
complicated than such a simple mean relation. 

In summary, the value of the optical absorption deter- 



mined by Greiner et al. ( 1995 ) is intriguing and signifi- 
cantly larger than what we obtained above: Ay = 3.49. 
On the other hand, their "confirmation" of the distance 
cannot be trusted, even if, mentioning the 3-5 kpc range 
of Tingay et al. ( 1995| l, they claim in their conclusions 



that: "The scattering of X-rays by the interstellar dust 
allows to derive a distance of GRO J 1655^40 of 3 kpc." 



4.2. Is there a systematic effect? 

The significantly larger value of A*~ Tdys compared to 
^optical j s j m p 0rtant smce ^ can provide strong insights 
into the role of dust, and the possible superabundance 
of dust over hydrogen close to the source itself. It is 
not only an important point for the distance determina- 
tion, but also for our understanding of the closeby en- 
vironments of such objects and therefore deserves to be 
verified. 



Jonker and Nelemans (2004i claim that optical ab 



sorption obtained from X-rays are systematically larger 
than those obtained from optical data, as shown in their 
Table 3, in which the values of A° ptlClil and A^" ldys are 
listed for 14 sources, among which GRO J 1655^40 
and A0620-00. Taken at face value, four targets in 
the list of 14 sources don't have a value in one of the 
two bands (an "X" is marked in its place): GS 1009- 
45, XTE Jl 118+480, H 1705-250, SAX J1819.3-2525. 
Moreover, five others have, strictly speaking, consistent 
values either because the values agree within the (large) 
uncertainties or because only upper/lower limits were 
determined (GRO J0422+32, A0620-00, GRO J 1655- 
40, GX 339^4 and GS 2 Section 023+338), leaving 
5 systems only for a true comparison: GS 1124-684, 
4U 1543^17, XTE J1550-564, XTE J1859+226 and 
GS 2000+25. 

4.2.7. GS 1124-684 

For GS 1124-684 (=GU Mus, Nova Mus 1991), 
leng et al. ( 1992) ) determine the color excess to be 
E(B - V) ~ 0.29, which corresponds to A° ptical = 0.9 ± 
0.1, by fitting a model of a HST/FOS spectrum be- 
tween 1600 and 4900A. 



The value of A 



X-rays _ 



1.28 



0.06 is computed by 
Jonker and Nelemans| (|2004|l taking N H = 2.28 x 10 21 



cm 



of 



Greiner et al. 



(1994a i and using the relation 



Nh/Av = 1.79 x 10 21 cm" 2 by |Predehl and Schmitt 



( 1995 [_ corresponding to a color excess of E(B -V) ~ 



0.4. But |Greiner et al.| ( |T99 4a ) write that their estimation 
of the color excess is confirmed by the works of, first, 
della Valle efaT1 ( [T99T] l, and second, |Chehg et al.|(^992> 



which is the reference for A 



optical 



Moreover, 



della Valle 



et al. ( 199 1) > use the problematic 6-extragalactic super- 
novae relation of Barbon ~e"t al.| ( |1990| l discussed above 
and adopt a mean value for the color excess derived 



7 It is actually misspelled "Schmit/f" in the caption of Table 3 of 
Jonker & Nelemans. 



from a range of values between 0.2 and 0.35, i.e. in 



agreement with that of |Cheng et al. ( 1992| l. 

As mentioned by Greiner et al. ( 1994a[ ) "the redden- 
ing of the X-ray transient in Muscae was first derived to 
E(B - V) « 0.2 - 0.3 from IUE and optical measure- 
ments", and give three references, among which two 
IAU Circulars by [West et aL] (|1991| but see also Ap- 



pendix [Aj, and Gonzalez-Riestra et al.| ( | 199 1| ). Unfor- 
tunately, none of the two circulars report a value of the 
color excess as such. The third reference given is a con- 
ference proceeding of Shrader C. and Gonzalez-Riestra 
R. 1991 (in "Workshop on Nova Muscae 1991", Lyngby 
1991, DRI Prep. 2-91, p. 85) which is unfortunately not 
referenced in the NASA ADS system. 

We can safely discard the source GS 1 124-684 in Ta- 
ble 3 of | Jonker and Nelemans] ( |2004| i. 

4.2.2. 4U 1543-47 

For 4U 1543-47 (=IL Lup*), the optical value A v = 
1.55 + 0.15 is quoted from |Orosz et al. ( 1998] > who also 
mention that their value is smaller than that derived from 
X-rays, quoting the values of Greiner et al. ( 1994b| l and 
|van der Woerd et al. (1989) which i s the r eference for 
A^" mys used bylJonker and Nelemans 



(2004). 



Strictly speaking, Orosz et al. (T998 see their Sec. 
4) have the following range 0.45 < E(B - V) < 0.55, 
while the values of the color excess obtained from X- 
ray data ranges from "0.56 to 0.77 assuming Ay — 
NhI 1.79 x 10 21 ( |Predehl and Schmittl|l995|l" That is, 
the two ranges almost agr ee. |Jonker and Nelemans 
(2004} ignore the value of Greiner et al. ( 1994b| l and 
choose the upper limit to derive A° ptlcal = 2.4+0. 1 . Given 
the difficulties of obtaining reliable estimates in both 



energy regimes, this disagreement between A 
^x-rays ^ ^ rea jjy conclusive. 



optical 
V 



and 



4.2.3. GS 2000+25 

For GS 2000+25 (=QZ Vul, Nova Vul 1988), the 



value from X-ray observations quoted by Jonker and 
Nelemans ( |2004) > is Ay = 6.4+1.0, citing |Tsunemi et al. 
( 1989). However, in this latter paper, the value quoted is 
A v = 4.41 (or \og m N H = 22.06 ± 0.006, see their Sec. 
EH, a, iii), said to be in agreement with th e optical esti- 

3.5) 



mation by Chevalier and Ilovaisky 



(1990 



• optical 



which is the reference given by Jonker & Nelemans for 
the optical value, who also note the large uncertainty of 
the latter. This source must therefore also be discarded 
from Table 3. 

4.2.4. XTE Jl 550-564 and XTE J 1859+226 

As for the two remaining sources, namely 
XTE J1550-564 (=V381 Nor) and XTE J1859+226 



10 



(=V406 Vul), there seems to be a difference between 



. X-rays 



and A 



optical 



but not necessarily because of a 
systematic effect in one of the two wavelength regime. 



For the X-ray binary XTE J 1550-564, Jonker and 



Nelemans 



(|2004| quote A° ptical = 



(1991) is pictured in their Fig. 6a where the non-zero 
slope showing the effect is caused by 2 deviating points 
only. According to their own caption, one point is a 
high-mass X-ray binary (4U 1516-56, a.k.a. Cir X-l) 



2.5 (no uncertainties) that we know today contains local absorption (Johnston 



computed from Sanchez-Fernandez et al. ( 1999) using 
their value of the equivalent width of the Diffuse Inter- 
stellar Ban d (EW=1.9A, f or the 4430A DIB) and the 
relatio n of Herbig (1975]). Sanchez-Fernandez et al. 
( 1999|) actually report A° p =2.2 from {E{B - V)> = 



0.7 ± 0.2 using a combination of the Nal lines and the 
DIB and using again the 6-supernovae relationship of 
Barbon et al.| ([T990) discussed above. Given the uncer- 
tainties on E(B - V), the two values actually agree. But 



if one uses the relationship of della Valle and Duerbeck 
( 1993}, E(B-V) = 1.6, and A° ptlcal > 5. In any case, their 



equivalent width of the Nal lines of 2.4A is much larger 
than th e significant range compu ted above from the 
study of jMunari and Zwitter| ( |l997} : EW Nal ~ 0.1-0.8A 
(or 0. < E(B - V) < 0.4). Moreover, they use a low- 
resolution (2A) spectrum, and it is not clear if the lines 



are saturated or not, which is the reason why Jonker and 



Nelemans (2004) use their value of the DIB only. 



As for XTE J 1859+226, [Jonker and Nelemans| l |2004) 
mention that A*~ rays =4.47 has no error bars. This value 



is quoted from the IAU Circular no. 7291 by dal Fiume 



et al. 

optical* 



dl999b who found N H to be "about 8x 10 21 cm ' 2 ". 



T.80: 



: 0.37 is given by Hynes et al. 

V) 



(2002 1 who 



found E(B - V) = 0.58 + 0.12 from a spectral fit of 
the UV feature at 2200A in their HST/STIS spectrum. 
Assuming A*' ldys being correct, the 2.6 mag difference 
is certainly due to local absorption in the source itself. 



et al. 200 1\ and where near infrared (NIR) spectra in 



the K band revealed itself to be too obscured to allow 
a spectral classification ( |Clark et al.| |2003| ) (see also 



Mignani et al. 2002 for HST data on this source). The 



second deviating point of the figure is a low-mass X-ray 
binary (4U 1728-337, a.k.a. GX 354+00) in a globu- 
lar cluster. The A° ptlcal of this source is determined by 



Grindlay and Hertz ( 1981| ) not with the mentioned opti- 



cal data but rather with JHK photometry. Moreover, the 
values of A^ ncd[ and A^" rays for this object agree within 
0.2 mag. 

We note that two sources (namely XTE J 1859+226 
and GRO J 1655— 40) seem to have A^" 1 ^ significantly 
larger than A° ptlcal . Given that most of the X-ray ob- 
servations are done during outbursts, and it is tempting 
to think that a local absorber makes the estimation of 



A*" rays flawed. 



We note however (hat the work of |Pre-| 



dehl and Schmitt ( 1995[ ), used by many for their relation 



between Ay and Nh, show precisely that there is a sur- 
prisingly good correlation between the simultaneously 
measured dust and gas density. 

It seems clear that this point should deserve more 
studies, especially with new methods of determining 



Xiang et al. (2007 ) use time delays 



^x-rays For instance, 
between photons that are being scattered by dust, and 
being recorded in the X-ray halo. Interestingly, these 
delays are also function of the distance to the source. 



4.3. Is there a discrepancy between Ay rays and A°y >lcal ? 

To summarize, the claimed systematic overestima- 
tion of the optical absorption obtained from X-ray ob- 
servations compared to the value obtained from opti- 
cal studies is not established, since only two sources 
(XTE J 1550-564 and XTE J 1859+226) among the 14 



used by Jonker and Nelemans (2004i potentially show 



any difference. To claim a systematic effect in one given 
direction also implies the knowledge of where the cor- 
rect value is. 

The overestimation of A*~ lays over Ay P " cal mentioned 



by Jonker and Nelemans ( 2004|) h as apparently already 
been observed by |Vrtilek et al.| (|1991|). But Vrtilek 



and coworkers obtained all their absorption values from 



optical data from an older source: Bradt and McClin 



tock (1983 1, who do not always measure the absorp- 
tion themselves but cite even older sources. Moreover, 



the supposed systematic overestimation of Vrtilek et al. 



4.4. The quiescence state might be variable 



In Jonker and Nelemans (2004), the quoted value of 



the absorption of GRO J 1655-40 is A 



X-rays 



= 4.8 ± 2.i 



w hich is still in agreem ent both with the peculiar value 
of 



Greiner et al. 



Hynes et al. 



(1998 



1995 i, and the value of A° ptIcal from 
A° p lcal =3.7±0.3) and our own deter- 



mination of A?? 103 =3.49 above. Nevertheless, one ad- 
ditional source of uncertainty is whether the object was 
observed in quiescence. Jonker and Nelemans ( 2004[ ) 
refer to |Kong et al.| ( |2002| ) concerning GRO J1655-40. 
These authors found a value of the X-ray luminosity 
measured with Chandra lower than that of Greiner et al. 
( |1995] t with ROSAT, and even lower than another mea- 
surement made between two large outbursts separated 
by one year, obtained with ASCA by |Ueda et al. ( 1998) 
and |Asai et al.| ( fT998[ ). Kong and coworkers note that 
their Chandra observation occurred about 4 years after 
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the last outburst, and may therefore measure the quies- 
cent X-ray level more accurately (although this is not 
necessarily true). In addit ion, [ Garcia et al.| ( |200T| ) em- 
phasize that Wagner et al. ( 1994| i observed variations of 
the quiescent X-ray flux in V404 Cyg by a factor of 10 
with ROSAT. 

Finally, the variability in the optical and near-infrared 
has been recently demonstrated by |Cantrell et al.| (2008 ) 
for when the object is considered in X-ray quiescence. 
They have studied optical and infrared photometry of 
A0620-00 spanning the years 1999-2007. They clearly 
demonstrate the existence of three different states (pas- 
sive, loop, active) in which only the passive one truly 
represent quiescence. Therefore, great care must be 
used when modelling optical lightcurves, since it will 
have a direct influence on the observed quiescent visual 
magnitude, and the determination of the Roche-lobe ef- 
fective radius. These issues are discussed in the next 
section. 

5. The characterisation of the orbit and the sec- 
ondary star of GRO J1655-40 

We have shown the various difficulties at determining 
a robust value of the color excess and/or the absorption 
value. We now turn to the other problematic term in 
equation [2] the determination of the absolute magni- 
tude of the secondary star, emphasizing that none of the 
studies mentioned in this section were aware of the re- 
quirement to isolate true optical quiescent lightcurves, 
as mentioned above. 

5.1. Comparison of height different studies of the 
orbital and secondary star parameters of 
GRO J1655^0. 
To obtain the absolute magnitude of the secondary 
star, we need to combine information on the tempera- 
ture of the star and information on its radius (see equa- 
tion [6]). There are various problems arising with this 
approach, which are all exemplified in Table [T] where 
we extracted the quantitative results on the orbital pa- 
rameters of GRO J 1655^40 and its secondary star ob- 



tained by van der Hooft et al. ( 1997 ), Orosz and Bailyn 



(fT9^7]l,[vander Hooft et al.|(|1998|l,|Phillips et al.|(|1999[), 
Shahbaz et all dl999|l, IGreene et al.l d2001b, IBeer andl 



Podsiadlowski 



( |2002] > and |Shahbaz| (f2003]>. All these 
studies use either optical spectra or CCD photometric 
lightcurves in the optical domain, or both. Some re- 
analyse data published previously in order to correct a 
specific assumption in the given analysis. 

In Table Q] we have extracted the relevant orbital and 
secondary star parameters, being very careful to choose 



the best quantitative results according to the authors. In 
the lower part of the table, we have, in addition, com- 
puted the projected rotational vsim (in km s _1 ) using 



the relationship of Home et al. ( 1986) when both the 



mass ratio q = m2/m\ (where m\ is the mass of the 
black hole; in Table [T] we use the inverse definition: 
Q = m\l 1112) and the secondary radial-velocity semi- 
amplitude K2 were availabl^] 

0.49<3 2/3 

vsinz = ^2(1 + q)- 



Q.6q 2 ^ +ln(l +4 1 / 3 ) 



(14) 



We emphasize that this relationship is valid when the 
hypothesis of synchronous rotation is verified. When 
the inclination angle was available, we computed the 
rotational velocity v rot . The Roche-lobe effective ra- 
dius (in R Q ) is calculated using the formula of Paczyhski 
( |1971| l (see also |Jonker and Ne lemans (2004) for a dis- 
cussion) from the period P and the secondary mass m^. 



Rl = /WheLobe = 0.234 P^ml' 3 



(15) 



with P the orbital period expressed in hours and m2 in 
solar masses. Given the Roche radius and the orbital 
period, we computed the equatorial velocity of a star 
with the same radius, which allow a cross-check of the 
synchronous rotation hypothesis. 

When a temperature was available, we also computed 
the so-called Roche luminosity, using equation [6j the 
Roche radius and assuming a solar temperature T = 
5800/f. We then calculated the absolute magnitude 
My from this Roche luminosity, using the absorption 
A y = 3.49 determined above for GRO J1655-40, and 
the apparent magnitude my from the considered data. 
Finally, we also computed the corresponding distance. 

5.2. Notes on individual studies. 

The paper by Bail yn et al.| ( |1995b| l is not discussed 
because it uses observations that are the first on this tar- 
get in the optical domain, and as such, provides only 
very few quantitative results. We also ignore the study 



by Hynes et al. ( 1998 1, who present a large dataset, 
and discuss extensively the problem of the extinction 
towards GRO J 1655^-0, but use basically all dynam- 
ical parameters from Orosz and Bailyn ( j!997[>. Other 
studies, such as Buxton and Vennes (2001 ) or Gonzalezl 



|Hernandez et al. (2008 ) provide updates to a small sub- 
set of the orbital parameters, but either have large uncer- 
tainties or refer heavily to values determined elsewhere. 



8 We note that |Shahbaz et al.Hl999| and |Shahbaz1j2003f also use 
this relationship but there is a typo in their formula: the 1/3 exponent 
is misplaced. 
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We think that none of the studies mentioned in this sec- 
tion obtained an accurate and complete parameter set 
about GRO J 1655^40. It is nevertheless insightful to 
excerpt specific issues that compromise the accuracy of 
the quantitative results. 



5.2.7. {van der Hooft et al.\^1997\ 

The main issue in |van der Hooft et aL] ( |1997| l is the 
fact that observations were taken during the outburst, 
and the overall magnitude increases by 0.9 mag dur- 
ing the course of the observations. Their folded /?-band 
lightcurve clearly shows two different levels of maxi- 
mum (see their Fig. 3). The authors finally determine 
the mass ratio by fitting their lightcurve with a basic 
eclipsing-disk model with X-ray heating and two es- 
timations of the X-ray luminosity communicated pri- 
vately by Harmon. Given the intrinsic variability seen 
in the data, it is clear that only estimations can be made. 
With the component masses and the orbital period, it is 
in principle possible to derive a value for K2 and cal- 
culate the rotational velocity. However, given the large 
uncertainties on the mass ratio and the secondary mass, 
it would be meaningless. We nonetheless note that, if 
the system is synchronous, the allowed range for the 
equatorial Roche velocity (70.5-86.2 km s~') is in dis- 
agreement with a projected (and hence minimum) rota- 
tional velocity of 93+3 km s^'determined spectroscop- 



ically by Israelian et al. ( 1999 1 and with which Foellmi 
|et al.| ( |2006| l totally agree. 



5.2.2. \Orosz and Balfyh]^997\ 

Orosz and Bailyn| ( |1997| l used optical spectra taken 



during the outburst (in 1995) and outside the outburst (in 
1996), as well as BVRI photometric lightcurves (only 
in 1996). Their spectroscopic data comprises the data 



published in Bailyn et al. ( 1995b[l. Despite the fact that 



the source was variable and the secondary star is fill- 
ing its Roche lobe, they fit a sine wave to their radial 
velocity data (see their Fig. 1). This issue is corrected 
later by Phillip s~et al.| ( |1999] but see below). In order 
to determine the inclination and the mass ratio, they fit 
their lightcurves with a model comprising 16 parame- 
ters, among which 9 are fixed, and where the spectrum 
of the star is approximated by a blackbody. This latter 
issue is corrected in Beer and Podsiadlowski (2002, also 
discussed below). 

Concerning the distance, if we take the observed 
value of the apparent magnitude V = 17.12, their lumi- 
nosity L = 46.6L , combined with E(B - V) = 1.3 and 
a standard R = 3.1, we obtain D = 3.1 kpc. However, 
if we take the Roche luminosity (hence assuming syn- 
chronicity, as did the authors) inferred from their pre- 



cise mass ni2, the orbital period P, and the visual ab- 
sorption value we determined earlier (Ay = 3.49), we 
obtain D = 3.8 + 0.2, i.e. larger than the upper limit de- 
termined by Hjellming and Rupen ( 1995] > with the radio 
jets. Moreover, we note that the value of the luminosity 
decreases by a factor of two in the corrected model of 
Beer and Podsiadlowski (2002, see below). 



5.2.3. {van der Hooft et al.\ $998) 

van der Hooft et al.| < fl998| ) attempted to fit their VRi 



lightcurves with an ellipsoidal model in which the sec- 
ondary mass ni2 and the inclination angle i are left free. 
However, they used a fixed mass function that has been 



determined with the early data taken by Bailyn et al. 
(1995b) during outburst. Once the mass function is 



known, the primary mass, and thus the mass ratio (re- 
quired for the model) are known. Moreover, they as- 
sume a bolometric luminosity (strangely called L opt in 
the article) L = 41 L Q taken from [Orosz and Bailyn" 



( 1997| l, which does not corresponds to the correct value: 
L = 46.6L . They also assume E(B - V) = 1.3 mag 
taken from the IAU Circular by |Horne et al.| ( fT996| l dis- 
cussed in Sec. 13.31 With such fixed values determined 
elsewhere from outburst data, it is unlikely that the out- 
come will be accurate. 

In Table 3, we note that the distance implied by their 
study on the distance cover the range 3.0^1.2 kpc, due 
to the fact that they have a very uncertain ma, hence 
mass ratio, combined with a poorly constrained inclina- 
tion angle. 



5.2.4. {Phillips etal]$199ty 

Phillips et al. ( 1999} reanalyse the outburst spectro- 



scopic data only of Orosz and Bailyn (1997), ignoring 
the photometric datasets. They show that X-ray heat- 
ing of the secondary surface can significantly modify 
the radial velocity curve. They finally obtain a semi- 
amplitude of K2 = 196 km s , i.e. much lower than 
the previous value, and therefore derive an "updated" 
value of the mass ratio. Using the the allowed range 



for inclination angle determined by|van der Hooft et al. 



(1998), they can compute a new mass for the black 
hole. This study certainly demonstrates the need to 
correctly take into account the X-ray irradiation when 
modelling the radial-velocity curves (or alternatively, 
re-emphasize the need to choose a truly quiescent op- 



tical lightcurve; see again Cantrell et al. 2008| ). But it 
does not provide accurate results for the GRO J 1655^-0 
system itself. 

We note in Table [T] that their results imply a smaller 
rotational velocity vsin/ = 75 ± 5 km s , which is 
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again in large disagreement with the measured value 
mentioned above. 

5.2.5. \Greene et al.\§200l) 



Greene et al. (2001 ) obtained photometric lightcurves 



in the optical (B, V and /, said to be "quantitatively in- 
distinguishable from that of Orosz and Bailyn ( 1997) 



and van der Hooft et al. ( 1998 1"). They also obtained 



for the first time near infrared data (in J and K bands) 
in order to have better constraints on the role of the ac- 
cretion disk. Moreover, their modelling is based on the 
code presented in Orosz and Hauschildt (2000) which is 
more sophisticated than what Orosz & Bailyn used. The 
authors mention that they make use of "all of the binary 
system observables and their uncertainties", mention- 



ing the 29 points of the radial velocity curve of Shahbaz 
et al. ( 1999 ) as their second dataset (the first one being 
their own photometry, and the third being the rotational 
velocity of Israelian et al.| (|1999|l). Table [T] illustrates 
clearly that the results of Greene et al. ( 2001) are in per- 
fect agreement with that of Shahbaz et al. (1999), but 
not with that of Beer and Podsiadlowski ( 2002 ) who use 
basically identical photometric data (see below). 

The reason of the discrepancy is provided by |Beer 
and Podsiadlowski ( 2002|l who discuss extensively the 



results of Greene et al. (2001 1. Similarly to the work 



of Orosz and Bailyn ( 1997 1, Greene and coworkers im- 
plicitly allowed for arbitrary offsets between the dif- 
ferent lightcurves, which are not fitted simultaneously 
(in addition to not using any extinction and distance 
information). Greene et al. concluded that the accre- 
tion disk does not contribute to the lightcurves. On 
the other hand, Beer and Podsiadlowski ( 2002 ) found a 
much poorer fit of the multi-color data when no disk is 
present, concluding that no fully self-consistent models 
can be obtained without a disk. 



5.2.6. \Shahbaz et al\(l999^ and \Shahbaz\$2003ty 

1999) analyse new optical spec- 



Shahbaz et al. 



tra taken during X-ray quiescence and find K2 - 
215 km s _1 . The data was then reanalyzed by Shah 



baz ( 2003 ) who developed a procedure to determine the 
spectroscopic mass ratio in interacting binaries only by 
modelling the secondary star spectrum only. 

There are various issues with this approach, but 
the main one is with the data. In the introduc- 
tion, the author emphasizes the numerous problems as- 
sociated with the developed procedure (intermediate- 
resolution spectroscopy only, shape of the rotational 
profile, wavelength-dependent limb-darkening, Roche- 
lobe filling and so on), but the solution proposed to 



reduce "these uncertainties" is to "determine the ex- 
act rotationally broadened spectrum from the secondary 
star in an interacting binary", i.e. make the best possi- 
ble model. It remains true however that even a perfect 
model cannot work well with medium quality data. 

As matter of fact, the spectral range is comprised in 
a very small region between 6320A and 6550A, where 
there are many iron lines, with a resolution of 4.2A. 
With an observed rotational velocity of 93 km s" 1 (see 
above), the line are heavily blended. However, an inter- 
mediate step of the procedure involves fitting the con- 
tinuum of the spectrum. This issue is crucial, since the 
procedure rely on the strength of the lines, which ob- 
viously depends critically on the true continuum level 
and the amount of possible veiling by the accretion 
disk. There is no detailed indication about this central 
step, apart that a spline fit is made before combining 
the spectra (which is also not an easy task). This issue 
is already present in the analysis of the abundances in 



GRO J 1655^-0 secondary star by Israelian et al. ( 1999 
and which has been questioned in Foellmi et al. ( 2007 



who demonstrate the impossibility to measure abun- 
dances as precisely as it has been claimecj^] Even with a 
Signal-to-Noise ratio of 100, it is very difficult to know 
where is the continuum in the spectrum, if any is truly 
present given the very broad lines. This requires a spec- 
trophotometric calibration of excellent quality (one or 
two % maximum), which is the case neither in Shahbaz 
|etal]([T999} nor in |Israelian etaLl ( [T999] l. 

Interestingly, Shahbaz ( |2003| l obtain a poor first fit 
that they attribute partially to the fact that some ele- 
ments might be overabundant as it has been supposedly 



shown by Israelian et al. ( 1999 1. Given the amount of 
uncertainties and the quality of the data used, we can 
safely discard the results of this study. Moreover, the 
distance implied by the results is the largest of the height 
studies compared in Table[T](D = 4.0 ± 0.5). 



5. 2. 7. Beer and Podsiadlowski ( 2002) 
|Beer and Podsiadlowski| (|2002|l reanalyze the pho- 



tometric data of Orosz and Bailyn ( 1997| l, and ignore 
their spectra. They developed a better lightcurve mod- 
elling by not assuming a perfect blackbody spectrum 
for the secondary star, which has strong consequences 
on the relative fluxes between the B and the V filter 
bandpasses. Moreover, they fit simultaneously in a self- 
consistent manner the lightcurves in the different pass- 
bands. One can see in Table [T] that their results imply 



'The issue has not been truly addressed in Gonza lez Hernandez| 
|et al.| (2008l, where it is still not demonstrated how the continuum 
location is determined. 



15 



a decrease of luminosity by a factor ~2, and a smaller 
mass of the secondary star by 40% compared to that of 
Orosz & Bailyn. 



The problem in the models of Beer and Podsiad 



|lowski| ( |2002| l is that the distance is said to be a free 
parameter, although it is not clear if it has been allowed 
to go as low as 1.0 kpc. In many places the authors 
claim that the distance of 3.2 kpc of Hjellm ing and Ru-"| 
|pen| ([T995 ) has been used "to tighten" their results, since 
the three main parameters (distance, color excess and 
temperature) are highly correlated, and can compensate 
for each other. But if D is basically fixed, what room 
is left for the other parameters? Interestingly, |Beer and| 
Podsiadlowski (2002 ) found it reasonable to use a dis 
tance of 3.2 kpc since the distance by [Hjellming and 



Rupen 



of 



Bai 



( 1995]> is consistent with the other determination 



yn et al.| (|1995a[>, |Greiner et aL] (|1995|l, |McKay 



and Kesteven| ( |1994| l and |Tingay et al.| ( |1995| l, although 



they also note that the existence of a significant differ- 
ence between the jet inclination (~ 85°) and the disk 
inclination angles makes the model of the wiggles in ra- 
dio jet data not necessarily appropriate. We have shown 
above how problematic are these references. 

Let us mention here that a distance of 1.0 kpc im- 
plies a velocity f3 = 0.28 and consequently an inclina- 
tion angle of the jets of ~ 71°. It means that at 1.0 kpc, 
GRO J 1655-40 would probably not have its disk mis- 
aligned with the jets to the contrary of what has been 
claimed (e.g. M accarone} |2002| l . 

Finally, we need to emphasize here that the first best- 



fit steady-state disk model in Beer and Podsiadlowski 



( 2002 1 is discarded by the authors because it was giv- 
ing a distance "much lower than 3.2 kpc". The work 



by Beer and Podsiadlowski (2002) is certainly the best 
multi-color modelling of the lightcurves of GRO J 1 655- 
40 to date. Unfortunately, it remains unknown what the 
solution of their model was with a much smaller dis- 
tance. 

5.3. Summary 

Despite the wealth of data, its variety and the even 
larger wealth of calculations made either on the data it- 
self, or on the interpretation of the quantitative results 
(evolutionary tracks, Monte Carlo simulations, x 1 dis- 
tance minimization between spectra, etc.), it is hard to 
determine which study actually provide accurate, and 
not only precise, results. It is partly due to the fact that 
each paper studies only a subset of all parameters, tak- 
ing the missing ones from other publications. 

To summarize, we have clearly demonstrated how 
difficult the modelling of an interacting binary such as 
GRO J 1655^-0 is, and that a complete and consistent 



modelling of the photometric and spectroscopic datasets 
taken during true optical quiescence is still lacking. 



6. A new estimation of the distance of GRO J1655- 
40 



Recently, Guver et al. ( 2008 1 have derived the dis- 
tance of the galactic neutron star 4U 1618-52 using the 
fact that red clump giant stars can be considered as stan- 



dard candles (see |L6pez-Corredoira et aL |2002| where 
the method is well explained). Following this method, 
we retrieved 2MASS JHK magnitudes of stars in a box 
of 0.45° x 0.45° around the position of GRO J1655^10. 
On the 14659 stars, we constructed a Color-Magnitude 
diagram (CMD) and selected the red clumps giants in 
a similar way to what is done in Guver et al. ( 2008 ), 
and as shown in Fig. [3] We obtained 6033 stars that 
we binned in magnitudes with bins of widths 0.5 mags. 
For each bin, we fitted a gaussian on the histogram of 
the number of stars as a function of the color J - K. 
Each gaussian fitted peak provides a value of the color 
J-K, which is then identified to the color of red clump 
giant Stars of the given magnitude bin (Fig. HJ. Identi- 
cally to |Guver et al. (2008 1, we then assume an abso- 
lute magnitude of K = -1.62 and a unabsorbed color 
(J - K)q = 0.7 for these stars. (See also |Durant and van] 
Kerkwijk 2006 , for slightly different values and a little 
variant to this method.) 

By using the general equation [2] the following re- 
lationship between the absorption magnitude in the K- 
band and the color indices: 



A K = c e x [( J -K)-(J- K) ] 



(16) 



with c e = 0.657 (see Guver et al.), we can derive a re- 
lationship between the distance and the absorption de- 
rived from the color, in this direction of the sky (see 
Fig.[5]l. We have also derived this relation with stars in- 
cluded in areas of 0.15° x 0.15° and 1.0° x 1.0° around 
the position of GRO J 1655-40(70527 stars in total and 
25108 stars selected for the box of one degree squared). 
It does not change the result significantly as shown in 
Fig. [5] We also made tests on the extent of the selected 
region in the CMD, and by varying the number, sizes 
and positions of the magnitude bins, but again it proved 
to have a negligible influence on the final curves. 
The problem is estimating the (J - K) and (J - K)q 



colours in equation 16 for GRO J 1655-40. The former 



can be visually estimated from the J and K lightcurves 
by|Greeneetal.| ( |200"T] l: J-K ~ 0.55 (|Beer and Podsiad- 



|lowski||2002| use J-K = 0.6). As for the later, we need 
to assume that the infrared color of the secondary star in 
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Figure 3: Color-Magnitude diagram of the stars in a square of 
0.45° X 0.45° around the position of GRO J1655-40 extracted from 
the 2MASS database. The two straight lines delimitate the selected 
region (red dots). The linked points indicate the maximum number 
density of selected stars in every magnitude bins. By using a larger 
area, we increase the number of stars and thus the contamination of 
the sample. On the other hand, using a smaller area make the results 
less precise. 
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Figure 5: Calculated relationships between the absorption Ak and the 
distance using the selected stars of the three different areas centered 
around GRO Jl 655— 40. The curves saturate around D ~ 5 kpc, or 
Ak ~ 0.4, certainly because of the increased amount of contami- 
nation of dwarfs stars and M-giants, as mentioned by Gu ver et al.| 
{2008} (non-completeness of the 2MASS data for faint magnitudes 
might also play a role, for the faintest ones). Similarly, at small dis- 
tance, there are too few stars to obtain a perfectly reliable curve. 
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Figure 4: Gaussian fits to the histograms of the number of selected 
stars for each magnitude bins as a function of the J - K color. One can 
see that there is a progression of the maximum of the gaussian peaks 
towards redder colors when we go from bright stars (lower curves with 
peaks to the left) to fainter stars (upper curves). 



GRO J 1655^-0 can be compared to that of an unaffected 
star of similar spectral type (F6IV). [Koornneef] ( |1983| l 
gave J - K ranging from 0.24 to 0.29 for F2 to F8 dwarf 
stars^l On the other, one can also use the 2MASS val- 



ues of the comparison star HD 156098 of Foellmi et al. 



(2006 ) which is 50 pc away from the sun and can be con- 
sidered as unaffected by extinction. Its infrared color is 
(J - K) = (J - K)o = 0.397. Using the extreme values 
of the given range 0.24 < (J - K)o < 0.4, and equa- 
tion 16 we derive a /T-band absorption 0.1 < Ak < 0.2 
for GRO J 1655^1-0. A comparison with the curves of 
Fig.|5]shows that the distance of GRO J1655^10 is cer- 
tainly less than 2.0 kpc. 

We note however that it is still not possible to def- 
initely conclude the exact distance of GRO J 1655-40 
(i.e. it lies at the same distance as the cluster NGC 6242) 
because of the impossibility to derive a reliable curve 
for the brightest (and thus nearby) red clump giant stars. 
However, these new results confirm that of Foellmi et 
al.: GRO J1655^0 is likely to be much closer than cur- 
rently admitted. 



No values are given for F giants, but J — K ranges from 0.22 to 
0.35 for F2 to F8 supergiants. 



17 



7. The distance of A0620-00 

We have presented above a confirmation that the dis- 
tance of GRO J 1655^-0 is certainly smaller than 2 kpc. 
One can thus ask: is GRO J 1655-40 the closest (stel- 
lar) known black hole to the Sun? According to Jonker 



Galactic plane (1=209.96°, b=6.54°; R.A.=6' ! 22 m 44.4 s , 
Dec.=-00°20 m 45 s ). 

More recently, |Gelino et al.| ( |2001| l determine a dis- 
tance of 1 164+1 14 pc for A0620-00 using the color ex- 
cess of |Wu et~aL (1976]). This distance has also been 



|and Nelem ans (2004), A0620-00 is at a distance of 
1.2±0.4 kpc from the Sun and is so far the closest 
known black hole to the Surp"| However, its distance 
is also problematic, but for different reasons than that of 
GRO J 1655^-0. We start by critically review the pub- 
lished distance of A0620-00 before presenting a new 
estimation of its maximum distance. 



Shahbaz et al. 



used by 

use the older value of 1.2±0.4 kpc citing 



( |2004l ). |Galloetal.|(|2006]) still 



Shahbaz et al. 



( [19941 ), |Gelino et aT] ( |2001| ) and Jonke r and Nelemans 
( |2004| ). Finally, we note that |Pal and Chakrabartr|(|2005[ ) 
still use the distance range of Shahbaz et al. ( 1994 1, and 
that|Esin et al.|(|2000| who obtained D - 1.4 kpc) cite 
|Shahbaz et al.|jl994| ) and |Barret et al.| ( [T996] ). All these 
studies, even those being very recent, directly rely on 
the 30-years old determination of E(B - V) by Wu et al. 



7.1. What is the the color excess towards A0620-00? even if the y do not cite the ori g inal P a P el " 



nes, 



In their study of distances of low-mass X-ray bina- 

quote 



Jonker and Nelemans 



(1994) and 
A0620-00. 



Barret et al. 



(2004) 
(1996) 



Shahbaz et al. 



for the distance of 



Shahbaz and collaborators give a distance 
range between 650 and 1450 pc, with a preferred value 
of 1050 pc. 

To estimate the extinction toward A0620-00 I Shah- 1 
|baz et alJ, ( |1994| > use an estimation of the color excess 
E(B - V) = 0.35 from |Wu et al] ( fT983] ). But Wu and 
collaborators actually quote their own results obtained 
a few years earlier: Wu et al. ( 1976| ). This latter pa- 



per describe that the extinction is measured by filling 
the "extinction spectral feature" at 2200A in their UV 
spectrum. However, the spectrum of Wu et al. (1976) 
has been obtained with the Astronomical Netherlands 
Satellite, and consists of no more than 5 points only, 
respectively at the central wavelengths of 1550, 1800, 
2200, 2500, and 3300A, simply because the instrument 
onboard the satellite had only 5 channels. One can ques- 
tion the reliability of such measurement given the ex- 
tremely low resolution, and the imperfect fit used to de- 
rive the color excess. 

Interestingly, |McClintock and Remillard| ([2000) pub- 
lished a HST/STIS spectrum of A0620-00 ranging from 
1900 to 3100 A, clearly revealing the absence of a "fea- 
ture" at 2200A (see their Fig. 1, upper panel). This ab- 
sence of any strong extinction feature in the UV spec- 
trum of |McClintock and Remillard| ( |2000| ) (who cite 
|Barret et al.[ |1996| for the distance) seems to show 
that the extinction might be particularly low. This is 
consistent with a location of A0620-00 far from the 



1 1 We obviously ignore here the possible problems with the deter- 
minations of all other microquasars listed in Jonker and Ne lemans] 
(2004) , and in particular XTE Jl 1 18+480 which has a published dis- 
tance of D= 1.8 ±0.6 kpc. 



7.2. A wealth of uncertain estimations 

|Barret et al.| ([T996) use two different methods to es- 
timate the distance of A0620-00, and found a value 
(1.2 kpc) "in agreement with previous determinations". 
We note that they do not apply either method to 
GRO J 1655-40 (which is also discussed in the paper) 
for which they take the literature value, 3.2 kpc. 

The first method is the magnitude-comparison 
method (equation [2]), assuming that the Roche-lobe ra- 
dius derived from dynamical studies is equal to the 
"radius" of the secondary star. Moreover, Barret and 
coworkers state that the absolute magnitude depends 
only on the spectral type. Using another method to 
compute the Roche radius requiring a mass assumption 
(M = 0.4M Q for mid-K secondaries), Barret and col- 
laborators confirm their distance value: D = 1.2 kpc. 
No uncertainties are provided, but they estimate them 



at about 25%. They also cite Haswell et al. ( 1993 1 and 
Shah baz" et al.| ( |1994| ), who provide the optical derre- 
dened magnitude and hence rely on the non-existent fea- 
ture measured by Wu et al. d!976| ). 

Jke and Greenstein ( 1977 ), cited by Jonker and Nele- 
s|( |2004| ) on their Table 3. discussed in Sec. [4] esti- 
mate the color excess using simultaneously the relations 
of|Spitzer|([T948l E(B - V) = 0.44 with a large uncer- 



tainty), Wampler (1966 


E(B - 


- V) ranging from 0.25 


to 0.6) and Aannestad and Purcell (1973). But Spitzer 


use the data of Stebbins et al. 


(|1940|), and Aannestad 


and coworkers use that oflYork 


(fT97T]).|Oke and Green- 



out much larger distances, since the object is largely out 
of the galactic plane". Finally, they use a distance of 
870 pc using another estimation of the absolute magni- 
tude of the secondary dwarf star, and claim that this new 
value is "somewhat smaller than previous estimates [...] 
but not inconsistent with them". 
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To summarize, it appears that all methods are more or 
less giving about the same result around 1 kpc, but with 
still a rather large scatter around this value. We also 
note that most estimations rely on the determination of 
the color excess by Wu et al. ( 1976| l, even in most recent 
literature. 

Below, we use VLT/UVES spectra to apply the 
maximal-distance method of Foellr m~et al.| (2006) 
which, even systematically uncertain, provides an 
model-independent estimation of the distance. Inter- 
estingly, we can measure the equivalent width of the 
sodium doublet in these UVES spectra of A0620-00: 
EW = 0.5 + 0.1 A. The value is slightly outside the al- 
lowed range given in Sec. |3.1| and more importantly, 



equation 1 1 might simply not be applicable to an object 
outside the galactic plane, in the direction of the anti- 
center, since it is outside the directions where the re- 
lationship has been calibrated. If we use it nonetheless 
with caution, we obtain E(B-V) = 0.22 + 0.07. Assum- 
ing R = 3.1, we have Ay — 0.68. According to Shahbaz 
et al. ( 1994 1, the secondary star radius is 0.8 R Q , its tem- 
perature 4000 K (implying a luminosity of L = 0.14L o , 
that is an absolute magnitude of My - 7.0), its visual 
magnitude V = 18.2 and 40% disc contamination (i.e. 
a magnitude increase of 0.55 mag). The corresponding 
distance is 1.3 kpc, and 1.6 without veiling. 

Could actually A0620-00 be further away? Given the 
various uncertainties, especially on the secondary star 
radius, this distance is certainly not conclusive. 

7.3. Astrometric estimations of the proper motion 

We looked in public archives for images that could 
provide an astrometrical estimation of the distance of 
A0620-00. Unfortunately, it has not been observed 
with an imaging technique by HST. We were finally 
able to find only acquisition images in the ESO archive 
with the following instruments: 3.6m/EFOSC2 (ESO, 
La Silla Observatory), VLT/ISAAC and VLT/FORS1 
(ESO, Cerro Paranal Observatory). Unfortunately, 
ISAAC frames proved to be useless because of the too 
small number of stars in the field of view, and were 
discarded. Finally, we also retrieved SuperCOSMOS 
images, which have a much poorer pixel scale but pro- 
vide the largest baseline in time. Table [2] summarize the 
properties of the three datasets. 

In order to check if the star has been ejected from 
a cluster with a runaway velocity, similarly to GRO 
J 1655-40, we performed astrometrical calculations with 
the above archives images. We chose the EFOSC2 im- 
age to be the reference since it had the sharpest PSF, 
and selected 40 bright isolated stars that were visible on 



all other frames. We computed the geometrical trans- 
formation map between the reference and the other im- 
ages using the fitted positions of these 40 stars. Finally, 
we computed the expected position of the target on the 
EFOSC2 image for a given transformation matrix and 
its actual position on this image. The difference in pix- 
els was transformed into arcsecond using the pixel scale 
of the EFOSC2 CCD. Finally, the results were divided 
by the time elapsed between the reference and the given 
image. 

It appeared that the best transformation fit was ob- 
tained between the FORS1 and EFOSC2 images, with 
an rms = 0.1508 and 0.1506 pixels in the X and Y di- 
rection respectively. This translates into an uncertainty 
of about 30 mas. The very small observed difference 
of the position of A0620-00 between the FORS1 and 
EFOSC2 images was clearly below this value: -10.4 and 
9.0 mas in the X and Y directions respectively. For the 
SuperCOSMOS images, the rms achieved in the trans- 
formation fit is 0.65 and 0.20 pixels, corresponding to 
a uncertainty of 436 and 134 mas respectively. The ob- 
served displacement was again inside the uncertainties 
of 71 and 136 mas. 

Taking the 2.03 years difference between the FORS1 
and EFOSC2 images, we obtain a maximum proper mo- 
tion of about 15-20 mas/yr, which roughly corresponds 
to that obtained with SuperCOSMOS images that are 
22.08 years apart from EFOSC2 ones. It is possible to 
look at the possible presence of a cluster of stars within 
a given radius, taking into account this upper limit of the 
astrometric motion of A0620-00 on the sky. Looking at 
a region of radius of about 3 degrees using SIMBAD, we 



found two clusters: [KPR2005] 22 (Kharchenko et al 



|2005) , and C 0619+023 ( |Collinder|[T93T) . While there 
is no information on the latter, the former has an esti- 
mated distance of 1.5 kpc. It is located 2.2808 degrees 
away from A0620-00, corresponding to a projected sep- 
aration on the sky of ~59.7 pc at 1.5 kpc. Given the up- 
per limit of the projected proper motion of A0620-00, 
one can estimate its minimal age to be about 367 000 
years. 

The current data does not allow us to derive strong 
constraints on the proper motion of A0620-00, and 
therefore on its possible origin. Combined with the 
systemic radial velocity (4 km s _1 according to 



Marsh 



et al. 1994[), it would be interesting to simulate possible 
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galactic orbits for A0620-00, similarly to what Mirabel 
|etaLl ( |2002] i have done for GRO J 1655^10. 



Table 2: Summary of the archival data used to derive astrometry of A0620-00. The instrument/project, observatory location, date of observations 
and bandpass (filters) are indicated as well as the Field-of-View in arcminutes, and the pixel scale (in arcseconds per pixel). The pixel scale of 
SuperCOSMOS plates are computed given a plate scale of 67.14 '/mm and 10/jm pixel size. The FOV of SuperCOSMOS images is originally 
6.4x6.4 degrees, but smaller regions can be retrieved electronically. 



Instrument 


Location 


Date 


Filter (name) 


FOV ('x') 


Scale ("/pixel) 


SuperCOSMOS/UK Schmidt 


Siding Spring Obs. 


16/02/1983 


I(RG715) 


up to 384x384 


0.6714 


VLT/FORS1 


VLT/Paranal 


7/01/2003 


V (V_BESS+35) 


6.8x6.8 


0.200 


La Silla/EFOSC2 


La Silla (3.6m) 


16-18/01/2005 


V (V#641) 


5.3x5.3 


0.31 (2x2 bin) 



8. The maximal-distance method of Foellmi et al. 
(2006) and its application to A0620-00 

8.1. Issues in the method 

In Foellmi et al. ( 2006| l a new method allowing to 
estimate a maximum distance has been presented and 
applied to the case of GRO J 1655^-0. It is based on the 
comparison between the calibrated spectroscopic fluxes 
of the secondary star and that of a companion star of 
similar spectral type and luminosity class which needs 
to be close enough to have a negligible absorption. The 
comparison is made with spectra obtained with the same 
instrument (VLT/UVES) configured with an identical 
setup. The main issue with this method is that one needs 
to make an hypothesis about the absolute magnitude dif- 
ference between the two stars. More precisely, if we call 
/ the ratio of the spectroscopic fluxes between the sec- 
ondary star in the microquasar, and the comparison star, 
|Foehmi et al.| (|2006 ) derive the following relation: 



5 log 



(DzJ_ 



+ M 2 - Mi > 



(17) 



where a is the (spectroscopic) absorption toward the tar- 
get (i.e. the average absorption within the narrow wave- 
length range given by the UVES spectra used to make 
the comparison). D\ and D2 are the distances of the 
microquasar and the reference star respectively. Obvi- 
ously, a must be null or positive. The difference be- 
tween the two absolute magnitudes Mi and M2 is how- 
ever impossible to evaluate. Or, in other wods, how dif- 
ferent is the absolute magnitude Mi compared to that 
of the reference star, given that the spectra are directly 
compared. For GRO J 1655-40, |Foellmi et aT] (p006 1 



allowed for a large range and made the hypothesis that 
Mi is comprised between M2 - 1 and M2 + 1. They 
finally found that D < 1.7 kpc by comparing flux of 
GRO J 1655^-0 to that of 4 nearby stars of similar spec- 
tral types. However, because the absolute magnitude is 
unknown, this method has a systematic intrinsic uncer- 
tainty. 

Since then, various authors mention this distance re- 
vision (Caballe ro-Garcia et al.| |2006[ |Takahashi et al.| 



2008 Joinet et al. 2008 ), but simply continue to use the 
canonical value. |Sala et al.| ( |2007| > used XMM-Newton 
data to derive a connection between the inner radius of 
the accretion disk and the distance. They argue that if 
the distance is less than 1.7 kpc and its mass is less than 
5 M Q , then the inner accretion disk radius will be inside 
the gravitational radius of the black hole. However, as 
noted by Combi et al. ( 2007 1, the mass and spin of the 
black hole are uncertain. If the black hole were rapidly 



rotating as suggested by McClintock et al. ( 2006 1 and 
its mass less than 5 M Q , the horizon of the Kerr black 
hole could then be well inside the inner accretion radius. 
Furthermore, |Sala et al.| ([2007 ) admit that an estimate of 
the inner accretion radius using XMM data is model de- 
pendent. 



Lasota ( 2008 ) mention that if we were to accept the 



new and smaller distance, the secondary in GRO J 1655- 
40 would not be filling its Roche-lobe (a similar remark 
is made in Caballe ro-Garcia et al.[ |2007| who do not 
provide however any reference about it) and in this case 
too we would have to invoke a mass transfer instability 
to explain the outbursts. But we have seen above that a 
complete modelling of the multi-color lightcurves with 
a good understanding of the disk emission and the ab- 
sorption is still lacking. 

8.2. Application to A0620-00 

In order to check the distance of A0620-00, we ap- 
plied two different methods. The first one is that of the 
red clump giant stars. Unfortunately, it proved to be 
useless because of the too small number of stars in the 
2MASS database in the direction of A0620-00, even 
taking a box of one degree around it. We were simply 
unable to identify a Red Clump Giant star branch in the 
CMD. It can certainly be explained by the direction of 
A0620-00 in the sky, which points away from the galac- 
tic center. The other method is that of Foellmi et al. 
(2006). To complement the study, we derive constraints 
on its proper motion using archival imaging data. 

We looked for VLT-UVES archival data on A0620- 
00, with the aim at comparing them to the UVES spectra 
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Figure 6: UVES spectrum of A0620-00 (lower and upper CCDs com- 
bined), smoothed with a boxcar of 1 1 pixels. Major spectral features 
are marked. The inset shows a close-up of the sodium doublet region. 
|Neilsen et al.|{2008) identify with caution the possible presence of 
HeI/15875 in this region. 



of the UVES POP databaseFksee Bagnulo et al. 2003 1. 



We have found 20 spectra of the target (ESO program 
ID 066.D-0157(A), P.I. Maeder), split in two different 
wavelength range: 4790 to 5755A, and 5840 to 6805A. 
The spectra were taken in 2000, December 5 (12 spec- 
tra), 17 (2 spectra) and 21 (6 spectra). The have been 
reduced and calibrated in exactly the same way as for 
GRO J 1655-40 (see |Foellmi et al] [20061 to which the 
reader is referred for a detailed description of the re- 
duction and flux-calibration of the spectra). The aver- 
age spectrum is shown in Fig. [6] Moreover, we also 
retrieved from the UVESPOP the spectra of three dif- 
ferent single stars with similar spectral types and lumi- 
nosity classes: HD 10361 (K5V), HD 100623 (K0V) 
andHD 209100 (K4.5V). 



As in Foellmi et al. (2006), we compare the flux- 
calibrated spectrum of A0620-00 with that of the com- 
parison stars (ignoring the Ha line, in emission in 
A0620-00 spectrum). The results are shown in Fig. [7] 
and summarized in Table [3] Given that a must be null 
or positive, we can see that the resulting maximum dis- 
tance is significantly smaller than 1 kpc, with a mean 
around 0.4 kpc. We emphasize here that we accounted 
for 40% veiling of the disk (i.e. we multiplied the / 
ratios in the table by a factor 0.6) before computing the 
distance. Any smaller fraction of the contamination will 
bring the object even closer to the Sun. 



http : //www . sc . eso . org/ santiago/uvespop/ 
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Figure 7: Relation between the spectroscopic absorption a and the 
distance, following the equation fTT] The condition that the absorption 
must be positive is represented by the horizontal line. The light gray 
area shows the systematic uncertainty of ± 1 magnitude, while the 
dark gray area shows an uncertainty of 10% on the spectroscopic flux. 
It shows that the likely distance of A0620-00 is probably smaller than 
0.5 kpc. 



9. Summary & Conclusions 

It has been clearly shown that despite a large num- 
ber of variants in distance methods, it is of prime im- 
portance to be extremely rigorous with assumptions 
and with the use of results made earlier and usually 
by others. The distance of microquasars, and in gen- 
eral of compact objects in our Galaxy, will benefit in 
a few years of the results of the satellite Gaia. Until 
then, we will certainly have to combine various types 
of information to infer good estimates to the distance 
of these objects. We mention here the work of |La-| 
zorenko et al. ( 2001) who have obtained, through as- 
trometric measurements using VLT/FORS1, a precision 
of 30 m/croarcseconds, similarly to what is expected 
for the VLTI instrument PRIMA, currently being com- 
missioned at the ESO Cerro Paranal Observatory. This 
method requires a rather simple observational setup, a 
large field, and a large number of stars in the field. 
These requirements are easily fulfilled in the case of mi- 
croquasars in the galactic plane like GRO J 1655^-0, al- 
though it might again be more difficult or impossible for 
A0620-00. 

From the present work, we can conclude that: 



• The use of the color excess measured using the 
sodium lines is risky. Not only one must have a 
spectroscopic resolution high enough, but its ap- 
plicability range is limited. Moreover, one must 
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Table 3: The three dwarf K stars used to compute the maximum distance of A0620-00 are summarized. The star's name, spectral type, distance 
(computed from the fflPPARCOS parallax), absolute magnitude, (B - V) color (obtained from SIMBAD) and the ratio of the flux-calibrated 
spectrum of A0620-00 with that of the star are indicated. The last column indicates the ranges of maximum distance D of A0620-00 obtained 
through the constraint of a > with the two limits of the absolute magnitude M\. See text for details. The quoted / flux ratios are multiplied by 
a factor 0.6 before estimating the distance to account for a contamination of the accretion disk of 40% according to Shahbaz et al. 1 1994|. The 
uncertainty on the distance values is ±0.2 kpc. The error on the absolute magnitude Mi is computed from the error on the HIPPARCOS distance. 



Star 


Sp. Type 


Distance 


M 2 


B-V 




/ 


Max. D\ (spec) 






(pc) 


(mag) 


(mag) 






(kpc) 


HD 10361 


K5V 


6.6+0.1 


6.7+0.2 


0.85 


(5.2 


+ 0.9) 10" 4 


0.24-0.59 


HD 100623 


K0V 


9.5+0.1 


6.1+0.1 


0.81 


(6.3 


+ 0.9) lO- 4 


0.31-0.77 


HD 209100 


K4.5V 


3.61+0.05 


6.9+0.2 


1.06 


(3.9 


+ 0.9) 10~ 4 


0.15-0.37 



be careful when choosing the relationship between 
the equivalent width and the color excess. 

There is no proven systematic overestimation of 
the optical absorption as determined from X-ray 
data compared to that inferred from optical data. 

The optical data also needs to be checked for being 
taken truly during optical quiescence or not. This 
is crucial for lightcurve models. 

The quality of the data and the uncertainties asso- 
ciated to it must be assessed rigorously. It is mean- 
ingless to apply a sophisticated model on data of 
limited quality. 

The maximum-distance method proposed by 
Foellmi et al. (2006) contains a systematic uncer- 



tainty, although its main hypothesis is sometimes 
implicitly used by others. 

About GRO J 1655-40: 

• We have shown that although the upper limit of 
3.5 kpc is a rather firm measurement, the value of 
3.2 kpc has never been really measured. 

• We have also shown that the lower limit of 3.0 kpc 
is questionable and actually based on a question- 
able assumption on the interpretation of absorp- 
tion lines in the radio spectrum. Moreover, we 
have challenged the relevance of comparing two 
radio spectra in a region where there HI clouds 
with anomalous velocities. 

• We have estimated that E(B - V) ~1.0 and A v ~ 
3.49. 

• The peculiar absorption determined from X-rays 
towards GRO J 1655^-0 possibly indicate the pres- 
ence of local dust close to the object. However, the 
X-ray flux in quiescence is also very variable. 



• The studies of the orbit and secondary star of 
GRO J 1655-40 rely very much on the quality of 
the model of the disk, and on the assumed fixed 
parameters, often extracted from other incomplete 
studies. As matter of fact, no true model of 
the lightcurves of GRO J 1655—40 combining the 
radial-velocity measurement and the multi-color 
lightcurves has been made while letting all the pa- 
rameters truly vary at the same time. 

• A new distance method using red clump giant stars 
applied to GRO J 1655—40 however confirms a dis- 
tance less than 2 kpc. 

• Although not proven, the new estimation of a 
smaller distance of GRO J 1655-40 strengthen the 
idea of its possible origin in NGC 6242. It would 
make GRO J 1655-40 one of the closest black holes 
to the Sun. 

• At 1.0 kpc, B = 0.28 and there would be no 
misalignment between the disk and the jets, since 
6»~71°. 

As for A0620-00: 

• We have shown that the published distance and 
many confirmations of it are based (not always ex- 
plicitely) on a single measurement of the color ex- 
cess made 30 years ago on a spectrum made of 5 
points, which is moreover in contradiction with an 
HST/STIS spectrum. 

• We have applied the maximum-distance method of 
|Foellmi et aLlpOOg] ) to A0620-00 and found that it 
could be indeed located much closer, to a distance 
of -0.4 kpc. 

• The example of A0620-00 illustrate that normal 
distance methods based on comparing the appar- 
ent and absolute magnitudes are difficult to apply 
to A0620-00 since the estimation of the extinction 
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cannot rely on relationship established in the galac- 
tic plane. For this reason also, the red clump gi- 
ant star distance method could not be applied to 
A0620-00. 

• We have also used archival images to infer an up- 
per limit of the proper motion of A0620-00, which 
appears quite small. Although we have found a 
cluster of stars 2.8 degrees away, the origin of 
A0620-00 remains uncertain. 

Finally, to the question "What is the closest black 
hole to the Sun?", our answer is the black hole X-ray 
binary A0620-00. 
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A. Bibliographical glitches 

We mention here the few bibliographic issues that we 
have encountered in studying the question of the dis- 
tance of GRO J 1655^10 and A0620-00. The Nature 



paper by Bailyn et al. ( 1995b) quote Arnett & Bowers 
in 1978 about the maximum neutron star mass, while 



the correct year is 1977 ( |Arnett and Bowers 19771. 
Phillips et al.| ( [T999] > quote "Bailyn et al. (1996b)" from 
which they use the outburst data. However the cor- 
rect year is 1995. Trimble and Leonard (1995) men- 
tion the IAU Circular nb. 6063 of |Inoue et al.| ( fl994"] l 



but attribute it for some reasons to Reynolds & Jauncy. 
When looking for the X-ray fluxes of GRO J1655-M 
in |Garcia et aT] ( [7998] ) quoted by|Mirabel et al.| (|2"002"li, 
the authors referred to |Barret et al. ( 1996| l. However, 
the value of the X-ray flux is actually quoted from an- 
other paper, said to be "in press", by Zhang et al. in 
1996, in the Astronomy & Astrophysics Supplement 
Series. There are three publications that could match 



this reference: Zhang et al. (1996a} that is describ- 
ing X-ray BATSE/CGRO observations through Earth oc- 
cultation and where GRO J 1655-40 is not mentioned, 
Zhang et al.|(|1996b|l that is about GRO J 1849-03 only 



and Zhang et al. " p996c[ l, in which only the source 
4U 1608-522 is discussed. Moreover, there is no pub- 
lications with Zhang as a first author in 1997 and in 



A&AS. There is actually a paper in 1997 in ApJ: Zhang 



eTaLl ( [1997) 1 where ASCA and BATSE X-ray fluxes of 



GRO J 1655^-0 are discussed, and this latter paper is 
certainly the correct reference. As for the color excess 
coming from other sources for GS 1 124-684 cited by 
|Greiner et al.| ( |1994a| l, there is a reference "West et al. 
(1991)". But in the reference list, we find only an IAU 
Circular with only "West R.M." as author, and West, 
R.M. Delia Valle M., Jarvis B., 1991, in "Workshop on 
Nova Muscae 1991", Lyngby, May 1991. Interestingly, 
the NASA ADS system does not list any of those two 
references. The only seemingly reference is the IAU 
Circular [West et al.| ( fl9"9Tj no 5165) but with the addi- 
tional and final coauthor Pizzaro G. not mentioned in 
the item of the reference list of Greine Fet al.| ([T994a). 
Finally, Orosz and Bailyn (19971 refer to a submitted 
paper by Robinson C. et al. (1996) in the Astrophysical 
Journal, that cannot actually be found in NASA ADS. 
We were actually not able to find any paper by C. Robin- 
son as a first, second or third author on this subject, in 
1995, 1996 or later. 
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